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Summary. — An attempt has been made to analyse some experimental 
data of high-energy showers so far reported. It is shown that the de- 
pendence of the meson multiplicity on the primary energy seems to be 
explained phenomenologically, at least within the present limitation 
of experimental accuracy, in terms of the single nucleon-nucleon collision 
by making allowance for the fluctuation in inelasticity. The result seems 
to be consistent with the assumption that the average energy of pro- 
duced mesons in the center-of-mass system is fairly constant (about 
200 MeV kinetic energy) independently of the primary energy at least 
up to about 100 GeV (in the laboratory system) so that the meson mul- 
tiplicity normalized at the value 1 of inelasticity increases rather sharply 
with the primary energy. At higher energies, the collimation of high- 
energy mesons into the forward and backward cones would result in 
increasing the average meson energy so that the increase of the normalized 
multiplicity with the primary energy seems to become slower than at 
lower energies. The results of the analysis are discussed in comparison 
with the existing theories of multiple meson production. 


1. — Introduction. 


There have been found by numerous workers the phenomena which could 
be explained in terms of single very high-energy nucleon-nucleon collisions 
in nuclear emulsions. Yet the plausible interpretation capable to apply to 
all the existing experiments is not established because of the experimental 
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inaccuracies and of the complexity of the phenomena. The dependence of 
meson multiplicity on the primary energy, for instance, could not yet be de- 
termined with a sufficient accuracy. As is well known, if we plotted the number 
of charged secondary particles against the estimated primary energy for these 
high-energy jets, the experimental points were distributed very widely. This 
fact has led some authors (1) to conclude that the observed jets consist of two 
groups with large and small inelasticities respectively. On the other hand, 
AMAI et al. (2) have pointed out recently the possibility that the apparent dis- 
persion of experimental points might be explained in terms of a possible dif- 
ference between the meson-nucleon and the nucleon-nucleon collisions. 

In order to examine the energy dependence of meson multiplicity, however, 
it seems probable that we must take into account the fluctuation in inelasticity 
of individual events. Accordingly, by confining ourselves to those events 
which could be well explained by the single nucleon-nucleon interaction and 
by making allowance for the inelasticity of the event, we may obtain a satis- 
factory interpretation of the energy dependence of meson multiplicity. In 
particular, an inference from the knowledge of recent experimental results. 
of the antiproton-nucleon annihilation process may make the situation at 
lower energies more plausible. With this aim we make an attempt to analyse 
some experimental data now existing. 


2. — Analysis. 


Although we are mainly interested in the single nucleon-nucleon interactions, 
we are forced to take those events which appear to be nucleon-nucleus col- 
lisions, because only a very few, well established examples of the elementary 
processes have been reported. Thus, we have selected the events satisfying 
all the following criteria: 


i) The number of heavily ionizing prongs emerged from the shower is 
less than four. — This criterion assures, at least at relatively low energies 
(less than several hundred GeV), that the interaction has taken place probably 
between the incident nucleon and a hydrogen nucleus or a single nucleon bound 
at the periphery of a heavier nucleus. At extremely high-energies, however, 
this criterion does not exclude the possibility that the interaction has occurred 
in a central region of the target nucleus, since the incident nucleon may drive 


(1) C. Ditwortn, 8. J. GoLDSACK, T. F. Hoane and L. Scarsi: Nuovo Cimento, 
10, 1261 (1953); T. F. Hoane: Journ. de Phys., 14, 395 (1953); G. BERTOLINO and 
D. PESCETTI: Nuovo Cimento, 12, 630 (1954). 

(?) S. Amat, T. Murota and M. NIsHIDA: Progr. Theor. Phys., 17, 807 (1957). 
C. DizworTx et al. (!), Hoane (!) and other authors have already suggested in a 
qualitative manner, the presence of this possibility. 
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a tunnel through the nucleus giving a very small excitation energy (*). But 
the present analysis is carried out tentatively on the assumption that this 
criterion ensures the event to be due mainly to a single collision between two 
nucleons. The event which bears some evidence for the existence of successive 
interactions of the primary particle inside the target nucleus, for instance a 
jet (3-+55,) observed by LOHRMANN (*), has been rejected. 


ii) At least for some of the secondary particles, the energies have been 
estimated on the basis of multiple scattering measurements or electron pair 
observations. 


ili) The Lorentz transformation from the laboratory to the center-of- 
mass systems has been carried out by taking into account the angles and 
energies of secondary particles. 


iv) The inelasticity of the interaction has been estimated from data ii) 
and iii). — The event in which the inelasticity was estimated simply from 
the maximum angle of shower particles has therefore been rejected. 


In Table I are shown the events (5-1) thus selected for the present analysis. 
It should be noted that the primary particles of these cosmic-ray showers were 
altogether singly charged except for shower 12. The number of heavily ionizing 
(black and grey) prongs accompanied by the shower is given in column 4, 
which does not involve tracks identified as slow 7-mesons. The multiplicity 
of charged mesons, n,, shown in column 6 involves these grey and black 
m-mesons, whereas it does not the identified secondary protons of relativistic 
velocities. For the event which did not allow to identify shower particles, we 


(3) F. C. RoesLER and C. B. A. McCusker: Nuovo Cimento, 10, 127 (1952); 
C. B. A. McCusker and F. C. RoESLER: Nuovo Cimento, 5, 1136 (1957). 

(4) E. LOHRMANN: Zeits. f. Naturfor., 11a, 561 (1956); Nuovo Cimento, 3, 822 (1956). 

(5) W. H. Barxas, R. W. BIRGE, W. W. CHuPP, A. G. Ekspone, G. GOLDHABER, 
S. GOLDHABER, H. H. HECKMAN, D. H. PERKINS, J. SANDWEISS, E. SEGRE, F. M. SMITH, 
D. H. STORK, L. vAN Rossum, E. AMALDI, G. BARONI, C. CASTAGNOLI, C. FRANZINETTI 
and A. MANFREDINI: Phys. Rev., 105, 1037 (1957). 

(6) S. KANEKO, 0. Kusumoro and S. Marsumoro: Journ. Phys. Soc. Jap., 11, 
1129 (1956). 


(7) E. LOHRMANN: Nuovo Cimento, 5, 1074 (1957). 

(8) K. GorrstEIN and M. TEUCHER: Zeits. f. Naturfor., 8a, 120 (1953). 

(@) L. v. LINDERN: Zeits. f. Naturfor., 11a, 340 (1956). 

(10) V. D. Hopper, S. Biswas and J. F. Dargy: Phys. Rev., 84, 457 (1951). 

(11) A. DEBENEDETTI, C. M. GARELLI, L. TALLONE and M. Vigone: Nuovo Cimento, 
4, 1142 (1956). 

(12) R. G. Grasser, D. M. HAsKIN, M. ScHEIN and J. J. Lorp: Phys. Rev., 99, 
1555 (1955). 

(13) P. Crok, M. Danysz, J. GIERULA, A. JURAK, M. Migsowicz, J. PERNEGR, 


J. VRANA and W. WOLTER: Nuovo Cimento, 6, 1409 (1957). 
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have assumed that two tracks among them were relativistic protons, Le., 
nN, =—2, where n, denotes the number of shower tracks. 

The data given by BARKAS et al. (5) are quite different in nature from the 
other listed in Table I. These authors have investigated extensively the anti- 
proton-nucleon annihilation process by the nuclear emulsion technique and 
have found that on the average (5.3 + 0.4) (charged and neutral) 7-mesons 
are produced in the process. The process is equivalent in available energy 
to the nucleon-nucleon interaction of the Lorentz factor y, = 2 and the in- 


elasticity Æ— 1 in the center-of-mass system (14), since these antiprotons came 


TABLE I. — Examples of high-energy interactions employed for the analysis. 


Average 
Sho- pecs Lorentz | Charged kinetic 
LA Primary heavily factor | meson | Inela- |energy of 
| Authors energy) |)... in the | multi- | sticity | mesons 
num-, ionizing os È i; 
neem (GeV) prongs, c.m.8. | plicity K in the 
| N Ve ny c.m.s. 
i (GeV) 
1 BARKAS et al. (>) — — 2.0) (4) 3.5 TOR) 0.18 
2 | KANEKO et al. (°) 25 0 SIT + 0.41 0.34 
3 | LOHRMANN (7) 37 0 4.5 17) 0 0.18 
4 GOTTSTEIN, TEUCHER (8) 40 0 4.7 18 1.0 083, 
5 LOHRMANN (°) 130 0 8.6 10 0.4 0.25 
6 LINDERN (°) 300 0 1183 8 05237) || Wear ©) 
7 LOHRMANN (?) 700 3 19 34 0.75 0.46 
8 HoPPER et al. (1) 1000 0 24 6 0.1 0.31 
9 LOHRMANN (7) 1700 0 30 10 0.31 1.0 
10 DEBENEDETTI et al. (1) 3800 2 45 38 0879 0.9 
11 GLASSER et al. (12) 20000 2 100° 13 0.12 PAL 
12 CIOK et al. (13) 330000 0 420 16 (4) 0.1 a few 


(a) The available energy in this antiproton-nucleon annihilation process is equal approximately 
to the mass energy of both of the particles; i.e., the process is equivalent to the case of the nucleon- 
nucleon collision of y.= 2 and K=1. 

(b) This event cannot be a pure proton-proton collision because n+ is not even. 

(c) These values have been estimated by us making use of the data of meson energies in the 
center-of-mass system given by Lindern. 

(d) These authors have regarded the jet 0 +14x as being due to a single nucleon-nucleon collision 
and have estimated the number of charged particles produced in the interaction to be 12. However, 
taking into account an anomalously large number of 7°-mesons as compared with charged mesons 


they have observed, we take the value 16 as a more probable estimate for the multiplicity of 
charged mesons. 


(4) This equivalence may be somewhat questionable because of the large inter- 
action cross-section observed for the former process. But the difference between them, 
if it exists, will be small so far as we are concerned with the process of meson pro- 
duction (see, for instance, N. Yasrma and K. KoBAYAKAWA: Progr. Theor. Phys 
19, 192 (1958)). 
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to rest or slowed down to very small velocities before they have annihilated. 
The results of this experiment provide us perhaps at present with the most 
reliable estimate of meson multiplicity in this low energy region because of 
their excellent statistical accuracy and of the well-defined energy. 

For the events shown in Table I we have plotted in Fig. 1 the normalized 
multiplicity n,/K, the multiplicity of charged mesons divided by the inelas- 
ticity, as a function of y, — 1 which represents a half of the maximum energy 
available for meson production in the center-of mass system in units of the 
nucleon rest energy. For the same data, on the other hand, the meson mul- 
tiplicity n, is plotted against the quantity K(y,—1) in Fig. 2. Of course, 
the latter quantity corresponds to a half of the center-of mass energy trans- 
ferred to meson field in the interaction. 


3. — Interpretation and discussion. 


3:1. — The values of primary energy and inelasticity estimated by each 
author for the cosmic-ray showers of Table I should be in considerable error 


Fig. 1. — The normalized multiplicity, the multiplicity divided by the inelasticity, of 
charged mesons plotted against the kinetic energy (in units of the nucleon rest energy) 
of the incident nucleon in the center-of-mass system for the events shown in Table I. 
The solid curve shows the straight line of slope 1 passing through the antiproton 
point (marked with a square). The broken curve is that calculated by the Heisenberg 
theory (Eq. (7)) assuming the $ power spectrum of produced 7-mesons. The dotted 
curve represents the straight line of slope 3 fitted to the higher energy points. 
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because of a severe limitation of the experimental technique and of uncer- 
tainties of the assumptions employed. Accordingly, we cannot draw the final 
conclusion only from a comparison between Fig. 1 and 2. It seems, however, 
that the experimental points shown in Fig. 1 can be fitted better to a smooth 
line not only than the usual plot of n, versus y, — 1 but also than the plot 
in Fig. 2. (A further discussion will be given later on.) 


100 nt 


Fig. 2. — The multiplicity of charged mesons plotted against the product of the inela- 

sticity and the kinetic energy (in units of the nucleon rest energy) of the incident 

nucleon in the center-of-mass system for the events shown in Table I. The figure 

attached to each experimental point refers to the shower number listed in column 1 

of Table I. The curve is the straight line of slope 1 passing through the anti- 
proton point (marked with a square). 


In particular, as shown in Fig. 1, the straight line passing through the 
antiproton point with slope 45° is fitted fairly well to the data at lower 
energies (less than a few hundred GeV). This line results evidently from the 
assumption that the average energy of produced mesons in the center-of-mass 
system is independent of the primary energy and is the same as that in the 
antiproton annihilation process. In this connection, it should be remembered 
that the antiproton data are in fact reliable. Furthermore, the cosmic-ray 
data in this energy region seem to be more accurate than those at higher 
energies, since the former events can perhaps definitely be regarded as due 
to pure proton-proton collisions except for shower 3 and also the energy measu- 
rements can be done with more ease in general at lower energies. Therefore, 
it seems reasonable to conclude, at least within the present experimental limi- 
tation, that the average energy of produced mesons is nearly constant inde- 
pendently of the primary energy at least up to about 100 GeV so that the 
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meson multiplicity shows a rather sharp increase with the primary energy; 
i.e., we obtain the empirical relation 


(1) Ne 223.0 K(y,—1), for ye< 7 


It is also shown that the experimental information obtained simply from 
the angular distribution of shower particles seems to support the above ex- 
pression. If we assume that the 7-mesons are emitted with a unique energy 
y. (less than y,) in the center-of- 
mass system, the maximum angle x D aerrenne'®) 
of shower particles in the labo- Ni + Hoang (17) 
ratory system, 0,,, is given ap- 
proximately by (1°) 


v—l\t_ Vr 
GER (G),;, = | ce ) Za 


(18) OY 


4 Pickup & Voyvodic we 


From a simple energetical con- 
sideration, we then obtain 


Fig. 3. — The quantity 9(y,— 1)/(y, sin 6,,) 
measured simply from the secondary angular 
The right hand side of this equa- distribution is plotted as a function of 


tion can be determined from the ye—1 for cosmic-ray showers observed by 
data of secondary angular distri- various authors. The straight line is the same 


bution by the usual method. (The as the solid curve in Fig. 1 (Eq. (1)). 
median angle method gives a sys- 

tematic error for estimating the value of y. in the case of y.>y,, but 
the error is very small (5) for y_ = 2.) In Fig. 3, we have plotted the 
experimental value of this quantity against y, — 1 for the showers observed 
by BERTOLINO (!), HOANG (!), and Pickup and Voyvonpic (1%), which in- 
volved showers initiated by nucleon-nucleus collisions (N, < 4). The fact 


(5) C. C. DirwortH, S. J. Gotpsack, T. F. Hoane and L. SCARSI: Nuovo Cimento, 
1 012 COSE 

(16) G. BERTOLINO: Nuovo Cimento, 3, 141 (1956). The showers with large inelas- 
ticities (belonging to zone A) observed by him have not been plotted here. Recent 
analysis made by I. A. IvanovsKaya and D. S. CHERNAVSKY (Nucl. Phys., 4, 29 (1957)) 
indicates that these showers were produced presumably by the collision of an incident 
nucleon with two or more target nucleons. 

(17) T. F. Hoane: Journ. Phys. Rad., 14, 395 (1953). 

(18) E. Pickup and L. Voyvopic: Phys. Rev., 84, 1190 (1951). 
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that the straight line of Eq. (1) fits fairly well to most of the lower energy points 
out of them implies 


Va =e sin bn = const (~ 2.4) ; 


as can be seen from Eq. (2). It should be remarked that these data which 
were obtained simply from the secondary angular distribution are in fair 
agreement with those which were analysed in more detail under the stringent 
selection criteria (1°). 


3-2. — At higher energies, deviations of the experimental points from the 
45° line can clearly be seen from Fig. 1 (Fig. 3 gives also a similar indication). 
This fact leads us to suggest that the increase of the average energy of mesons 
resulting from increasing primary energy takes place at about 100 GeV. The 
dependence of the average energy upon the primary energy suggested by this 
view is as a matter of course in substantial agreement with a general trend 
of the values estimated directly by the respective authors which are shown 
in column 8 of Table I. 

This is also consistent with experimental results on the angular distribution 
of produced mesons. Many of the experiments so far performed by means of 
photographic emulsions and cloud chambers show a nearly isotropic distri- 
bution up to about 100 GeV primary energy (**8:92°) and some degree of aniso- 
tropy at greater energies (21). At extremely high-energies, a strong anisotropy 
has been well established (118). It seems, therefore, that the increase of the 
average energy is closely connected with the appearance of the angular ani- 
sotropy. 

This view would also be reasonable from theoretical considerations. A 
possible strong interaction of produced x-mesons with each other may result 
in a fairly constant value of the average energy of a number of mesons with 
an isotropic angular distribution at lower primary energies, whereas at higher 
energies (y. — 7) the Lorentz contraction of the interaction volume may play 
an essential role resulting in a collimation of high-energy mesons into the 


(1°) By analising his own data, BERTOLINO (!) has obtained the empirical relation 
n,— 1 = 3.33K(y,— 1) + 0.975, 


which is quite similar to Eq. (1). 

(2°) For example, R. Graccont, A. LovatI, A. Mura and C. Succi: Nuovo Cimento, 
4, 826 (1956); Y. Warase, S. Mirani, K. Sua and Y. TANAKA: Progr. Theor. Phys., 
18, 314 (1957). 

(21) S. ToxunaGa: Journ. Phys. Soc. Jap., 12, 454 (1957); W. WINKLER: Helv. 
Phys. Acta, 29, 267 (1956). 
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forward and backward cones, as was explained by HEISENBERG (222%) or 
LANDAU (#4), in addition to a nearly isotropic ejection of lower energy mesons. 
In fact, this expectation is consistent with many experimental evidence (712) 
that a large fraction of secondary particles have been emitted with relatively 
low energies even in energetic collisions. 


3°3. — In order to compare the present results with the theories of multiple 
meson production, let us confine ourselves only to the data at higher energies. 
It appears that the straight line of slope 4 is fitted fairly well to showers 5 
to 12 as shown by the dotted curve in Fig. 1; i.e., we have empirically 


(3) Mi Ci ye ahi ’ for Ve To 
with the tentative values of 


029 
(3a) 
me 0.5. 


On the other hand, any straight line in Fig. 2 is written 
(4) n, = O'[K(y.—1)]”, 


where O' and m' are constants. If we assume that only m-mesons are pro- 
duced in the interaction, it is evident that Eq. (3) holds in the case when the 
average value of the total energies of the 7-mesons in the center-of-mass system 
is given by 


(5) ie ce (92 = ies ; 


i.e., when it depends primarily on the Lorentz factor of the colliding system. 
On the other hand, if the average energy depends on the energy transferred 
to the meson field in such a manner as 


(6) Vr oc [K(y. a 1) 67 1 


then the multiplicity is given by Eq. (4). It is possible that the main dif- 
ference between Eqs. (3) and (4) can be given by the physical contents of 
Eqs. (5) and (6), even if other particles (heavy mesons and nucleon pairs) 


(22) W. HEISENBERG: Zeits. Phys., 126, 569 (1949). 
(23) W. HEISENBERG: Zeits. Phys., 133, 65 (1952); Kosmische Strahlung (Berlin, 


1953), p. 148. 
(24) S. Z. BeLEN'KyJI and L. D. LANDAU: Suppl. Nuovo Cimento, 1, 15 (1956). 
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having different average energies from that of x-mesons come to play a role, 
since a contribution of these particles to the multiplicity is likely to be fairly 
small even at extremely high-energies (?°). 

As already mentioned, the present data do not permit us to judge defi- 
nitely this alternative, nevertheless we intend to give a preference to Eq. (5) 
by the following reasons. First, from Fig. 2 an inference could be drawn that 
the higher energy points (say, showers 9 to 12) are considerably below the 
straight line of slope 45° which is well fitted to the lower energy points 
(showers 1 to 5), so that a smooth line capable to be fitted to all the experi- 
mental points in the whole energy region could never be drawn. Furthermore, 
the assumption implied by Eq. (6) that the average energy is determined by 
the transferred energy alone rather than the Lorentz contraction of inter- 
acting fields would be less plausible on theoretical grounds. Lastly, an other 
support for this preference might be found in the fact that several examples 
of meson showers of very high multiplicity have been observed thus far at 
relatively low energies (426), since these events might otherwise not be ex- 
plained even if we assumed the occurrence of cascade processes inside the 
target nucleus. It seems probable, therefore, that in the higher energy region 
the average energy of mesons in the center-of-mass system is determined 
primarily by the Lorentz factor of the colliding particles so that the energy 
dependence of multiplicity is given approximately by Eq. (3). 


4. — Comparison with the- theories of Heisenberg and Landau. 


On the point of view described above, we have to call attention to the 
Heisenberg theory of multiple meson production (2223). This is because in his 
theory the maximum and hence the average energies of a particular kind of 
mesons are determined by the Lorentz factor of the incident particle and also 
because the inelasticity is introduced explicitly into the formulation. The 
theory can be characterized by the power law energy spectrum of mesons 
produced: 

dn = dy, /yn'*. 


If we assume only the production of z-mesons with «= 2, for example, the 
multiplicity of charged mesons can be written approximately (22) 


(7) m= 4.5K(y.—1)1—y,*)(yi—1)>. 


(25) F. A. BRISBOUT, C. DAHANAYAKE, A. ENGLER, Y. FUJIMOTO and D. H. PERKINS: 
Phil. Mag., 1, 605 (1956). 


(7°) M. ScHEIN, D. M. Haskin and R. G. Grasser: Suppl. Nuovo Cimento, 12, 
355 (1954). 
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This equation is shown by the broken line in Fig. 1, which is not inconsistent 
with the experimental results. We may obtain à better fit by taking into 
account the production of heavy mesons and nucleon pairs. However, the 
inverse square law (x = 1) of the energy spectrum (23) seems to give too high 
a multiplicity at the highest energy. 

A comparison between the present results and the Landau theory (21) can 
be made under the assumption that the completely inelastic (K — 1) inter- 
action corresponds to the case of a head-on collision in the theory. The theory 
is apparently in good agreement with the fact that the empirical value of the 
exponent m is approximately } (Eqs. (3) and (3a)), but the numerical co- 
efficient C~1 predicted by the theory is too small in comparison with the 
empirical value C~9. However, by introducing the effective inelasticity as 
evaluated by AMAT et al. (27) into the theory, we may obtain a better agreement 
with the experiments. Therefore, the present analysis could not permit us 
to discriminate between the theories of Heisenberg and of Landau. 


5. — Conclusion. 


The results obtained above are summarized as follows: 


i) By making use of the normalized multiplicity of mesons produced, 
we may explain consistently various high-energy showers in terms of the single 
nucleon-nucleon interactions. 


ii) At lower primary energies (from 10 to about 100 GeV), it seems that 
the average energy of the produced mesons in the center-of-mass system is ap- 
proximately the same as that in the antiproton-nucleon annihilation process 
(about 200 MeV kinetic energy) independently of the primary energy and the 
inelasticity. Thus, the normalized multiplicity increases sharply with the 
primary energy (Eq. (1)). 


iii) At energies greater than 100 GeV, say, the average meson energy 
seems to increase with the primary energy resulting in a slower increase of 
multiplicity. These circumstances; have, at least to some extent, relation to 
the experimental evidence that the angular distribution of mesons is nearly 
isotropic up to about 100 GeV primary energy and it also shows a strong 
anisotropy at extremely high-energies. 


iv) Though any decisive conclusions should be reserved because of 
experimental uncertainties, the energy dependence of multiplicity at higher 


(27) S. Amar, H. Fuxupa, C. Iso and M. Sato: Progr. Theor. Phys., 17, 241 (1957). 
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energies may be approximated empirically by Eqs. (3) and (3a). The Heisen- 


berg and the Landau theories, in a somewhat modified form, are shown to 


be not inconsistent with the present results. 
* Kk OK 


We wish to express our sincere thanks to Professor Y. WATASE for his 
kind interest and encouragement. We are also indebted to Prof. S. HAYAKAWA, 
Prof. K. Suga, Mr. O. Kusumoto and Mr. 8. MATSUMOTO for their stimulating 
discussion and criticism. 


RIASSUNTO (*) 


Si analizzano alcuni dati sperimentali su sciami di alta energia finora resi noti. 
Si dimostra che, tenuto conto delle fluttuazioni dell’anelasticità, appare possibile spie- 
gare fenomenologicamente la dipendenza della molteplicità mesonica dall’energia pri- 
maria, almeno nei limiti dell’esattezza attuale della sperimentazione, in termini della 
singola collisione nucleone-nucleone. Il risultato sembra accordarsi con l'ipotesi che 
l'energia media dei mesoni prodotti è abbastanza costante (energia cinetica circa 
200 MeV) indipendentemente dall’energia primaria, almeno fino a circa 100 GeV (nel 
sistema del laboratorio), cosicchè la molteplicità mesonica normalizzata al valore 1 
dell’anelasticità cresce piuttosto bruscamente con l’energia primaria. Ad energie supe- 
riori, la collimazione dei mesoni di alta energia nei coni anteriore e posteriore risulte- 
rebbe nell’accrescere l’energia mesonica media, cosicchè l'aumento della molteplicità 
normalizzata con l'energia primaria sembra esser più lento che alle basse energie. Si 
confrontano i risultati dell'analisi con quelli delle attuali teorie della produzione mul- 
tipla dei mesoni. 


(*) Traduzione a cura della Redazione. 
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Photons, Gravitons and the Cosmological Constant. 


A. PERES 


Israel Institute of Technology - Haifa, Israel 


(ricevuto il 28 Novembre 1957) 


Summary. — An attempt is made to solve Maxwell equations in a Rie- 
mannian space subjected to Ry, = Agy: It is shown that the equation 
of motion for an electromagnetic wave-packet is similar to that of a 
particle (photon) having rest mass (f,/A/c). We then deal with the 
equation R,, = Ag,, itself, and seek for a solution having the form of 
gravitational waves. It is found that, if such waves can carry quanta 
of energy, then the rest mass of an elementary wave packet (a graviton) 
is (4y/—2A/c). It may therefore be concluded that A cannot be negative, 
and if gravitons do exist, it can also not be positive. It is therefore 
necessarily zero, and the rest masses of photons and gravitons are also null. 


1. — Electromagnetic waves in curved space ('). 


We seek a solution of Maxwell equations in vacuo: 


tk im kn 
=9 9 Fan 


Pn = Aia pa Aia 
the space-time metrics being subjected to 


ka = AGix ; 


(1) A less rigorous derivation of the same result has been given by the author in 
Compt. Rend., 239, 1023 (1954). 
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i.e. in the approximation of weak electromagnetic fields. We have 


RG =) = 


or, lowering index 1: 


But (2) 


= g A ik T 9" hik;i an APR) =e AN LA A,E; 5 


k 
A si;k Ta sas an 
In the last expression, the first term is null; because of Lorentz gauge con- 


dition ae. = 0, and the second is 


A,R° = AA 65 = AA,. 


Therefore 


This vector equation is very difficult to solve directly, because we cannot. 
set a system of cartesian co-ordinates, and derivatives of the /7 field will 
appear on the left hand member. Moreover, in curved co-ordinates, as the 
vector A, moves away, it undergoes an apparent change of direction which 
prevents us from defining a velocity for its components. As a matter of fact, 
the apparent velocity can be different for the various component of A,, and 
it can also be modified by infinitesimal changes of the co-ordinates. 

It is therefore better to seek for the equation of motion of some character- 
istic scalar tied to A;: The simplest such scalar is of course the magnitude. 
AN Oi 2158 


ea 
We have 
giA;A,; 
Ud 
id A A / k 
47 è "4 gd (7) 


(2) Much care is to be taken of the signs of the curvature tensor components, as 
they vary according to different authors. In the present paper, we shall follow ToLman 
Relativity, Thermodynamies and Cosmology (Oxford, 1934) and define 


i LE î i 8 i s a 
Bin = Lise — Lien + Din To Din Ta 
and 
i 
Fo = Rix: 


Ji 
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But 
DÉS Se es Ag” A,A; =— AA’, 
thus 


Ù x A 
Ag + AA = gt A; i : 


We can now choose, at the point which we are considering, a system of 
almost cartesian co-ordinates, where the /7 ocally vanish; they will also be 
vanishingly small in its meighborhood: as a matter of fact, if we consider a 
wave packet the dimensions of which need not exceed a few wavelengths, 
the /’ terms will have order of magnitude AA in the region of interest. 

Let us then seek a solution A,= B;f, where the B; are slowly varying 
functions (*) and f is a scalar (presumably satisfying some wave-like equation). 


We have 


4 
A 5% aa By, "la B;3] = Bite TR jk) 
A; RI B; Bin PÉDALE: lira 
1) =G),= B B3 i 
where B= \/ g"" By Bn. It is easily seen that 
B;\7 
g'iB; (7) == 0 , 
so that 
ste ii 5 Bix B; 
47 + AA = Ag ge? DA ER 


Thus, if we can find slowly varying 5, (*) the right hand side of the last 
equation will have order of magnitude A/?A?, i.e. will be vanishingly small 
even compared with the AA term in the left-hand side. That such B; can 
be found may be shown by developing the vector wave equation A,.,}"+ AA4;=0, 
in which we replace A, by A(B,/B). We then have 


D: B,\* B)\ * BD 
Ties 24 | ) 4( ) A ie =; 
i (7) NTS ae (ee B i 


(*) Le. in B;,, = B;; — B,T;; both terms on the right-hand side are of the same 
order of magnitude BAA. 
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or, with the help of the just derived value of 4,” 


B. Bay B, sk i TAI (7 au Hi ik 
DI 1j 2; k 9 = — 0 
ae B (7) eee a ARE. 


Each term of this equation involves covariant derivatives of B;, so that 
none is of zero order in AA if B, indeed has the form B;=C;+D,;, where 
the 0, are arbitrary constants and the D, are small terms of order C;22A. There 
remain only first (and higher) order terms, and the four equations can be solved 
for the four D; in a way consistent with the hypothesis that they are small. 

Once this point has been settled, we can go on with the scalar wave- 


equation 


Ar iA 0. 


This equation is known to represent waves travelling with group-velocity 


where » is the frequency. 
And if we have such a wave packet containing only one photon of energy 
W= hy we can consider its rest mass as being 


ce? Q70e È 


su | (IA hVA 


Up e? = 


2. — Gravitational waves in curved space. 
We shall now deal with the equation R,,= Ag, itself. We have 


OF ex ee 
dx" dar 2 dai 


1 i @ 
e 5 GE (Gee Tee) D (9""L'.mn) + terms quadratic in J’. 


As previously, we choose geodetic co-ordinates, such that the Z' vanish at 
the point which we are considering, and have magnitude AA in a neighborhood 
of dimensions Z We can therefore neglect terms quadratic in 7. 

Moreover, we can impose on the field four additional conditions (‘) such as 


Or ee AIX 


(4) Cfr. EinsTRIN: The Meaning of Relativity (Princeton, 1953), p. 87. 
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(Note that these conditions are to hold exactly, and not only approximately, 
in finite region of space-time. This is quite similar to the use of Lorentz 
gauge condition A4*,= 0 in the previous section). It then remains 


n o? dir 
al Sam age I 


As usual, we write gx =yx+h%x, where y; = +1, 0 are the galilean values 

for gi; and A, are small quantities, the square of which can be neglected. 
Then 

mn ir 

27 Bam ba 


= A(yix + hix) . 


Let us now see how a perturbation $h,, of the h;, propagates. (This per- 
turbation has to be small, of course, relative to the y,, but not necessarily 
with respect to the h,,;). 

We have 


Il mn = œ(h ik fa Sh;x) 


2? dae da” 


= Alix + hin + dhix) . 


We can substract from this the equation for the unperturbed %,, (which 
represents the « natural » curvature of the field) and get 


02(Shix) 


Oam 0x" 


= Al (Ohix) . 


1 pa 
2 


Once more we get a Klein-Gordon-like equation, which corresponds to 
waves travelling with group velocity 


2Ac? 


? 
472? 


pee 


and if there do exist wuch wave-packets carrying quanta of energy (gravitons) 
W = hy, their rest mass will be (5) 


= LA o hV—2A 
LEE CAR 2 


(6) This result had already been found, in a less direct way, by TONNELAT (cfr. 
DE BROGLIE: Théorie générale des particules à spin (Paris, 1954), p. 191). 

Much care should be taken of signs: our £,, is equal to the R,, there, but our g,, 
is minus the g,, there. Therefore the definitions for the A’s are opposite in signs, 
which explains the difference in the final formulae. 


35 - Il Nuovo Cimento. 


1129 


538 A. PERES 


3. — Conclusions. 


As it is a fact of experience that electromagnetic waves can carry quanta 
of energy, we may conclude that À cannot be negative (unless we bring some 
modifications to Maxwell or Einstein’s equations). If it is not null, its value 
is less than 5-10-51 em-? (5), so that the rest mass of the photon is certainly 
less than 10-* g. = 

On the other hand, it has still not been proved that gravitational waves. 
can carry quanta of energy. If they can, then A cannot also be positive. It is 
therefore necessarily zero, and the rest masses of photons and gravitons are 
also null. 


* Kk OK 


The author is indebted to Prof. N. Rosen for valuable suggestions and 
clarifying discussions of the results. 


(6) Cfr. TOLMAN: Relativity, Thermodynamics and Cosmology (Oxford, 1934), p. 474. 


RAS SUN ON) 


Cerchiamo una soluzione delle equazioni di Maxwell in uno spazio Riemanniano 
per il caso R,,= Ag;,. Dimostriamo che l’equazione del moto di un pacchetto d’onde 
elettromagnetiche è simile a quella di una particella (fotone) con massa a riposo #/4/e. 
Trattiamo poi l'equazione R,,= 4g; e cerchiamo una soluzione in forma d’onde gravi- 
tazionali. Si trova che, se tali onde possono trasportare quanti di energia, la massa a. 
riposo di un pacchetto d’onde elementare (gravitone) è 44/—2//c. Si può pertanto con- 
cludere che A non può essere negativo e che, se i gravitoni esistono, non può neanche 
essere positivo. E pertanto necessariamente zero e le masse a riposo dei fotoni e dei 
gravitoni sono a loro volta nulle. 


(*) Traduzione a cura della Redazione. 
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Cosmic-Ray Evidence on the Origin of Meteorites (*). 


S. F. SINGER 


University of Maryland - College Park, Maryland 


(ricevuto il 19 Gennaio 1958) 


Summary. — The birth of a meteorite coincides with the breakup of its 
parent body and also with the beginning of the meteorite’s exposure to 
cosmic rays. This continued exposure creates through nuclear evapor- 
ations definitely calculable amounts of *He, about 3 of it by way of 
tritium. Hence the date of breakup can be obtained from the *He con- 
tent; or better still, from the *He and tritium content measured in the 
same sample. However, the latter method has only limited applicability. 
The two methods give good accord; of six meteorites dated, all are found 
to be about 300 million years old. This agreement establishes confidence 
in the calculations of the *He production rate and leads to some con- 
siderations concerning the process of nuclear excitation. The age is also 
in good accord with Opik’s calculation of the mean lifetime of meteorite- 
like bodies against capture by the inner planets. From this agreement 
conclusions can be drawn concerning the mechanism of meteorite creation 
and concerning the relative constancy of cosmic-ray intensity as far 
back as 300 million years. 


1. — Introduction. 


It seems fairly well established, especially from the metallurgical evidence 
on iron meteorites, that meteorites were formed by slow cooling and under 
high pressure (12). Suitable conditions prevail only in the interior of good- 
sized bodies (planetoids) which may have existed at one time between the 
orbits of Mars and Jupiter in the present asteroid belt. According to current 


(*) Presented as an invited paper at the 1957 Spring meeting of the American 
Physical Society; supported in part by the Office of Scientific Research, U. S. Air Force. 

(1) H. H. UHnLIG: Geochim. et Cosmochim. Acta, 6, 228 (1954): 

(2) E. J. OPIK: Irish Astron. Journ., 3, 206 (1955). 
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theories, mainly associated with the names of KuIPER and UREY, these planet- 
oids were formed at the beginning of the solar system, first as dust aggregates; 
they condensed under gravitational attraction, were heated by radioactivity, 
and liquified. Within a very short time, by geological standards, these small 
bodies cooled. In this process the stony components rose to the top while the 
iron-nickel mass remained at the center. This differentiation, which is also 
thought to exist in the earth, gives rise to the two main types of meteorites, 
stones and irons. In the slow cooling of the core, the characteristic crystal 
structures of iron meteorites were formed showing the so-called Widmann- 
stàtten figures. Some time after, these planetoids were broken up by col- 
lisions and produced what we now call the asteroids. The term « meteorites » 
is applied to those asteroids which have fallen on the earth or on the moon. 
Experimental evidence, both from the lead isotope method (*) and from the 
argon-potassium method (‘) of radioactive dating, is now accumulating, showing 
that the solidification of the planetoids took place around 4-19° years ago; 
this corresponds to about the age of the solar system. 

A separate problem and one to which the present paper is devoted is the 
date of breakup of the planetoids, i.e., the date of creation of the meteorites. 
The method we use depends on the fact that up until breakup the (potential) 
meteorites are in the interior of the planetoid and are shielded from cosmic 
rays; but after breakup they become exposed to cosmic rays. The cosmic-ray 
exposure time thus gives us a method of dating the creation of a meteorite (5). 


2. — Methods of dating. 


21. Helium-3 Dating. — After it was realized that cosmic-ray bombardment 
of meteorites would produce substantial and measurable amounts of *He, 
the possibility of using a meteorite as an integrating cosmic-ray meter became 
immediately apparent (5). Dividing the calculated production rate of *He 
into the measured amount of *He, one immediately obtains the time interval > 
during which the meteorite has been exposed to the cosmic-ray flux (assumed 
to be of constant intensity), and therefore the date of creation of the meteorite 
from the breakup of its parent planet. This approach has been used to give 


(*) C. PATTERSON, H. Brown, G. TILTON and M. INGHRAM: Phys. Rev., 92, 
1234 (1953). 

(4) E. K. GERLING and T. G. PavLova: Dokl. Akad. Nauk USSR, 77, 85 (1951); 
G. J. WASSERBURG and R. J. HAYDEN: Phys. Rev., 97, 86 (1955); S. J. THOMSON and 
K. I. MAYNE: Geochim. et Cosmochim. Acta, 7, 169 (1955); W. GENTNER and W. KLEy: 
Zeits. Naturfor., 10a, 832 (1955). 

(5) S. F. SINGER: Astrophys. Journ., 119, 291 (1954). 

(6) S. EF. SINGER: Nature, 170, 728 (1952). 
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so-called « cosmic-ray ages » for a number of meteorites (7), or at least a lower 
limit. It was found for example, that meteorites with the highest *He content 
had breakup ages of about 300 million years, or possibly more (Table I). 


TABLE I. — Compilation of age determinations for various meteorites. 
Preatmo- ‘ Duration of 
Total He spheric 3He content | Estimated cosmic-ray 
SONORI diameter «| (107° em?/g) | « geometric exposure 
(1076 cm8/g) (em) factor » (yr) 
Mount Ayliff ana 14.6 8.80 DES 200 -- 300 : 105 
Tamarugal 23.6 42 5.57 4 SUI = 340-108 
Carbo 22.0 48 4.5% = 1 ~ 300-106 
Toluca (Durham) 18.9 Very large Ab Byes © nl =, 290-108 
Thunda 29.0 Small but 6.8? 1-2 200 — 400: 108 
irregular | 
Norton County [ 2.27 +0.11 ° 240-106 © 
(Stone) ] sam 
Furnas County | fall 2.35+0.11 280-106 © 
(Stone) [ 


(a) Reference (*). 
(b) P. REASBECK and K. I. MAYNE: Nature, 176, 733 (1955). 
(c) Reference (1). 


But to maintain this view, it had first to be shown that the low solidification 
ages (of about 100-200 million years) derived from the *He-uranium method (8) 
were invalid, either (i) because of the absence of uranium (°°) or (ii) because 
of the leakage of radiogenic ‘He (711.12), 


(7) S. F. SINGER: Sci. Amer., 191, 36 (1954). 

(8) J. DALTON, F. A. PANETH, P. REASBECK, S. J. THomson and K. I. MAYNE: 
Nature, 172, 1168 (1953). 

OMAN CG] Urny: Nature, 175° 321 (1955): 

(19) G. W. ReED and A. TURKEVICH: Nature, 176, 794 (1955). 

(1!) S. F. SINGER: Phys. Rev., 105, 765 (1957). 

(12) A crucial experiment (!!) has been suggested to test hypothesis (i). If it 
turns out to be correct, then we suggest the following interpretation of the U-+He data: 
the meteorite fragments which are.formed in the break-up are heated throughout by 
the shock wave of the collision. In many meteorites this heating is not high enough 
to destroy the crystal structure; but it is sufficient to accelerate the diffusion of the 
radiogenic helium which is located in the boundaries of the crystal grains, and thereby 
to give an apparent recent solidification age. This suggestion appears to us less ad hoc 
than the one put forward by P. REASBECK and K. I. MAYNE (Nature, 176, 186 (1955)), 
namely that the iron phase solidified only about 200 million years ago, and that the 
breakup occurred soon after this solidification. 
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The calculation of the *He production rate in a particular meteorite, and 
therefore the derivation of its breakup time, is made difficult by the fact that 
we do not know the shape and size of the meteorite before it enters the atmo- 
sphere, nor do we generally know the position of the sample with respect to 
the preatmospheric boundary. In most cases, therefore (e.g., for the meteo- 
rites Tamarugal and Toluca in Table I), we can only give a lower limit to the 
age. However, for the meteorite Carbo, helium determinations have been 
made as a function of depth (1), and for this meteorite we can estimate the 
geometric factor (14), and therefore the age, quite well. 

In many cases we can only say that the sample came from the interior of 
the meteorite and therefore experiences a production rate less than that at 
the maximum, i.e., we estimate a geometric factor of less than one. For very 
small meteorites, however, the cosmic-ray flux could have penetrated over a 
larger solid angle; here we must estimate a geometric factor greater than one. 
For example, for the meteorite Mt. Ayliff we estimate a geometric factor of 
the order of 2 to 3. The uncertainties, while large, do not change our conclusions 
substantially. 

The uncertainties arise as follows: 


1) It is very probable that Mt. Ayuirr did not enter as a small spherical me- 
teorite as for example PANETH (8) calculated. It may have formed part of the skin 
of a larger meteorite from which it was spalled off. In that case the geometric factor 
would be only one. This seems to have happened in the case of the Toluca meteorite: 
the sample at Durham has a very high *He content while.the sample in Hamburg 
has a ‘He content less than 1% of Durham (8), clearly indicating that the Durham 
sample came from the region of the skin whereas the Hamburg sample came from the 
interior. We therefore feel justified in assuming that the geometrie factor for Toluca 
(Durham) is not very much less than one, and that the age cannot be very much 
greater than 300 million years. 


2) But there exist in addition a large number of meteorites whose He content 
is very small, in some cases unmeasurable. We cannot at the present time tell whether 
these form the interior portions of larger meteorites, i.e. a geometric factor much less 
than one, or whether they are meteorites which were created by breakup only a rela- 
tively short time ago. In other words, we do not know very much yet concerning the 
dispersion of ages and the results of Table I should not be used to indicate a unique 
age for all meteorites. 


3) A further uncertainty of course, is the accuracy of the calculated value (5) 
for the *He production rate. Other investigators have derived much lower rates than 
ours and find therefore correspondingly higher breakup ages. For Mt. Ayliff FIREMAN (15) 


(5) F. A. PANETH, P. REASBECK and K. I. Mayne: Nature, 172, 200 (1953). 

(4) Defined as the ratio of *He production rate in the sample to the rate at the 
transition maximum, i.e., (2--15)cm from the preatmospheric boundary of a very 
large iron meteorite. 

(15) E. L. Fireman: Phys. Rev., 97, 1303 (1955). 
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quotes an age of 4000 million years while CuRRIE, Lippy and WoLFGANG (1) give 
its age as 1400 million years. In both cases the He production rate was calculated 
from experiments on the production of tritium in eyclotron-irradiated targets. (Pos- 
sible reasons for the discrepancy are discussed below). 


22. Tritium-Helium-3 Dating. — In a recent experiment BEGEMANN, GEISS 
and Hess (17) have applied the following method: in recently fallen stone 
meteorites (Norton County and Furnas County) where a measurable tritium 
activity existed, both the tritium and *He contents were measured. From 
the tritium content they can calculate the production rate for tritium, and, 
by derivation, for *He (Tables II and III). These measured production rates 
are in good accord with our calculations, as discussed below. Again under 
the assumption of constant cosmic-ray intensity, they derive a cosmic-ray 
age of about 260 million years, a value which is in general accord with our age 
determinations (7) (see Table I). 

It is evident that the tritium-*He method overcomes the major disadvantage 
of the *He method, namely the uncertainty of the position of the sample rela- 
tive to the preatmospheric boundary. On the other hand, there are very few 
cases of recent falls to which the tritium-*He method can be applied. It seems 
to us, therefore, that the determination of *He as a function of depth (1%) can 
be used most widely as a means of «locating » the sample and refining the age 
estimates. Ideally, the depth-variation method should be checked by means 
of the tritium method and a consistent calibration obtained. 


3. — Discussion. 


The good accord between the age deduced from the *He-tritium method 
and from our *He calculations is of course very gratifying. One needs to 
explain, however, (1) why the direct measurements (1516) of tritium production 
in a cyclotron are so much lower, and (2) how one one can explain the fact 
that the meteorite ages are as low as 300 million years. (3) Finally, we wish 
to make some remarks concerning the prehistoric cosmic-ray intensity. 


8 1. Rate of production of He and Tritium by cosmic rays. — In our initial 
publication (*) we calculated the rate of production of He (as given in Table II), 
and of tritium (!8) (as given in Table III) by the following method: 

The meteorite was considered to be of infinite size and the production rate 
was determined at the position of the transition maximum, about 2 to 15 cm 
deep for an iron meteorite. A quantitative calculation was carried out using 


. A. CURRIE, W. F. Lippy and R. L. WoLFrGaNG: Phys. Rev., 101, 1557 (1956). 
. BEGEMANN, J. Gæiss and D. C. Hess: Phys. Rev., 107, 540 (1957). 
. F. SINGER: Phys. Rev., 90, 168 (1953). 
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TABLE II. — Helium-3 production rate in meteorites. 


Production rate 
cms (STP) 
per g per year 


Location 
(relative to 
preatmospheric skin) 


Method 


to transition maximum 


(a) Reference (5). 
(b) Reference (!). 


eso KO At transition maximum Calculated from cosmic-ray data | 
and evaporation theory © 
~ 0.9-10-14 Unknown, but probably close | Calculated from measured tri- 


tium activity ? 


TaBLe III. — Tritium production rates in meteorites. 
| ductic t Activity : 
et e ì 3 Location Method 
| (atoms/g/s) (counts/min/g) 

Selo 0.50 At maximum near skin | Calculated from cosmic- 
of very large meteor- ray data and evapor- 
ite. ation theory ©. 

| STORE Ve 0.02 Reduces to large mete- | Calculated from expe- 
orite case by dividing riments using 2 GeV 
by a geometrie factor protons ?. 
of 4. 

Sisko ORs 0.05 At maximum near skin | Calculated from expe- 
of a very large me- riments with 430 MeV 
teorite. and 2 GeV protons. 

4.7-10-3 0.28 Unknown. Measured 9. 

(a) Reference (1). 

(6) Reference (15). 

(c) Reference (1°). 


(d) Reference (37). 


the present cosmic-ray flux near the earth, taking into account the transition 
effect due to low-energy cosmic-ray secondaries, and including also the inter- 
action of z-mesons which would not decay but would produce additional nuclear 
disintegrations. This calculation was checked by a more empirical approach 
in which we analyzed a large sample of nuclear evaporations, about 12000 stars 
obtained by the Bristol group, as to their excitation energy. Then, by ap- 
plying the results of the evaporation theory (1°), the rate of production of *He 


(19) K. J. Le CouTEUR: Proc. Phys. Soc. London, A 63, 259 (1950). 
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and tritium was determined. For the distribution of excitation energies under 
consideration about 65% of the production of *He is by way of tritium, while 
the rest is directly produced. More recently, a rate of production of tritium 
has been calculated based on experiments in accelerators with 430 MeV and 
2 GeV protons (1516), These values have been consistently lower than ours (18) 
(see Table III) without any obvious reason for the discrepancy. We cannot 
say for certain that the results of BEGEMANN et al. (1?) check our calculations; 
we do not know the position of their sample with respect to the meteorite 
boundary. Furthermore, there is some difference between stone and iron in 
the types of nuclear evaporations produced. However, their experimental 
determination of the cosmic-ray production rate for tritium seems to favor 
our higher value; in any case their value should be less than the value at the 
transition maximum. The diserepaney between our value and the value de- 
rived from the cyclotron experiments may perhaps be due to the fact that 
the excitation energy of the evaporating nuclei is much higher in the cosmic-ray 
case than for the cyclotron experiments because of the higher mean energy 
of cosmic rays. A 2 GeV proton will tend to give its energy to a knock-on 
nucleon and furnish only about 50 MeV for nuclear excitation (2°). Part- 
icles of the order of 10 GeV, however, would create many mesons in a single 
interaction, but in a heavy nucleus a substantial fraction of these mesons are 
immediately absorbed. In each reabsorption approximately 50 MeV will be 
given to the nucleus. The process of meson production and reabsorption within 
the same nucleus is a means of making available a large amount of energy for 
nuclear excitation. This hypothesis, if true, would mean that the cross-section 
for tritium and *He production should rise beyond 2 GeV. It would therefore 
be very important to continue the experimental work beyond this energy. 


3'2. Distribution of meteorites ages. - We can now inquire why the meteorite 
breakup times should show a clustering near 300 million years. Even though 
the sample is very sparse (6 meteorites) and the uncertainty in age rather high, 
the probability of such a close coincidence is very low and points to an under- 
lying reason. 

One can, of course, assume that all meteorites were created 300 million years 
ago but this seems hardly probable. A more likely explanation is one based 
on the work of ÒPIK who has shown (2!) that asteroidal bodies crossing the 
earth’s orbit have a mean life of the order of 100 million years (Table IV). 
Therefore, the probability that meteorites could have survived for much more 
than 300 million years is rather small. 

Starting with a few planetoids the rate of creation of meteorites probably 
first rose exponentially, as the numbers of fragments increased and therefore 


(20) I am indebted to Professor G. PuPPI for discussion on this point. 
(21) E. J. OpiK: Proc. Roy. Irish Acad.. A 54, 165 (1951); Contr. Armagh Obs., no. 6. 
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TABLE IV. — Lifetimes of objects having multiple crossing with inner planets ©. 


: Apollo Adonis Hermes Icarus jee Geminids 
Object P Encke's (meteor 
| Ce asteroids ==) comet shower) 
| | | 
| Inclination 6.4° 1.6° | 479 | 23.2° 12.5° 23.5° 
Eccentricity 0.57 Osiris | 0.47 0.83 0.85 0.90 
Mean life, 105yr | 72 72 | 42 | 170 270 250 
| (a) From reference (?!). 


also the probability of collisions. At the same time, the disrupting collisions 
also placed the fragments into eccentric orbits so that they were rapidly swept 
up by the planets. This constant attrition of meteorites by capture probably 
slowed the rising curve and let it reach a maximum. From then on, the number 
of new meteorites created by fragmentation must have dropped radically. 
At the present time, fragmentations of large bodies probably occur only occa- 
sionally, and can no longer be described by a continuous curve (on a time 
scale of 100 million years) but rather as random events. In this way the exis- 
tence of a peak in the vicinity of 300 million years can be interpreted quite 
reasonably as a particularly effective creation of meteorites in the last 500 mil- 
lion years. Earlier creations would have been attenuated so severely that 
they are no longer incident on the earth. 

We can say very little about the dispersion of the meteorite ages. Because 
of the reasons given earlier, we do not know how to find unambiguously the 
age of meteorites having a very small *He content. They could be regarded 
as having been created in much more recent times, say in the last 50 million 
years, without contradicting the present experimental evidence. The number 
of meteorites reported (5) to have a low *He content consistent with such an 
age is approximately 10, while the number having no detectable helium content 
is about 5. 


3°3. Remarks on the prehistoric cosmic-ray intensity. — We can make use 
of the good accord between Opik’s calculations of lifetime and the meteorite 
ages determined from the cosmic-ray method to make some definite remarks 
about the prehistoric cosmic-ray intensity. It gives one confidence that the 
cosmic-ray intensity has not varied greatly in prehistoric times. While we 
cannot, of course, exclude fluctuations of any magnitude, the average value 
probably was within a factor of 2 of its present value. 

The analysis of MARTIN (??) according to which the cosmic-ray intensity 
must have been 3 or more times its present value, is questionable (1!) since 


(22) G. R. MARTIN: Geochim. et Cosmochim. Acta, 3, 288 (1953). 
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= 


it is based on a solidification age for the planetoids of 100 million years obtained 
with the use of the ‘He-uranium method (823). 

A more precise statement would be possible if we knew the exact orbital 
characteristics of the meteorites under consideration. Using Opik’s theory (21), 
it would then be possible to calculate the probable lifetime of the object against 
capture by the inner planets. Assuming that the orbit resulted when the 
meteorite was created, namely that the collision which creates the meteorite 
also makes the orbit eccentric enough to intersect the earth, we would then 
obtain the precise exposure time to cosmic rays. 

The use, in addition to tritium, of other radioactive nuclides, e.g., 5*Mn 
with a suggested half-life of ~ 10° years (21), would help to check the con- 
stancy of cosmic-ray flux. This fascinating aspect of meteorite and cosmic-ray 
science is just beginning. 


(23) Neither can one conclude that the cosmic ray intensity has been 1/10 to 
1/20 its present value just because the *H :*He age is so much less than the 4A : 40K 
age (17). The two «ages » refer to different physical events. 

(24) R. K. SHELINE and J. E. HooPER: Nature, 179, 85 (1957). 


Note added: 


Our attention has been called to a paper just published by E. L. FIREMAN and 
D. SCHWARZER (Geochim. Cosmochim. Acta, 11, 252 (1957)). Their paper was submitted 
in August 1956, and represents therefore the first successful application of the *H-*He 
method to the measurement of the cosmic ray age of meteorites. The authors have 
examined two recently fallen iron meteorites and obtain the following ages: Para de 
Minas (1700-108 y.), and Norfolk (900-106 y.); both are much higher than the ages shown 
in Table I. This difference, however, only serves to reinforce our hypothesis (see point 2 
of Discussion) that planetary breakups occurred frequently throughout the history of 
the solar system. We would expect that the observed distribution of meteorite cosmic 
ray ages is the product of two functions: (1) the distribution of break-up dates (showing 
a peak early in the history of the solar system); and (ii) the exponential sweep-up 
factor (with mean life (100-200): 106 y.). 


RIASSUNTO (*) 


La nascita di un meteorite coincide con la frantumazione del corpo che gli dà ori- 
gine nonchè con l’inizio della sua esposizione ai raggi cosmici. Tale esposizione conti- 
nuata crea attraverso evaporazioni nucleari quantitativi esattamente calcolabili di *He, 
circa $ dei quali passando per il tritio. La data di nascita del meteorite può quindi 


(*) Traduzione a cura della Redazione. 
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essere desunta dal contenuto di *He; o meglio ancora dal contenuto di He e di tritio 
misurati nello stesso campione. Tuttavia, quest’ultimo metodo è di applicabilità limi- 
tata. I due metodi danno risultati in buon accordo fra di loro; sei meteoriti esaminati 
furono tutti trovati vecchi di circa 300 milioni di anni. Tale accordo convalida i cal- 
coli del tasso di produzione dell’ 8He e porta a qualche considerazione sul processo 
d’eccitazione nucleare. L'età riscontrata è anche in buon accordo coi calcoli di Opik 
per la vita media dei corpi meteoritici rispetto alla cattura da parte dei pianeti interni. 
Da questo accordo si possono trarre conclusioni sul meccanismo di creazione dei meteo- 
riti e sulla costanza relativa dell’intensità della radiazione cosmica negli ultimi 200 mi- 
lioni di anni. 
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Associated Production of Y° and K,.. 


M. S. SINHA and C. R. MALLIK 
Bose Institute - Caleutta 


(ricevuto il 5 Febbraio 1958) 


Summary. — A direct evidence for the associated production of a X° 
along with a Kj, has been obtained. The most interesting feature of 
the picture is that, besides showing the decays of a A%+n-+p) and 
Kj(>-y*+») it shows the materialization of the photon resulting from 
20 A°+y. The values of the energy of the emitted photon as derived 
from the decay dynamics and also from the energies of the soft electrons 
produced by it agree remarkably well. An overall energy balance is 
obtained for the whole event which is represented by the following reaction 


Tiara 


Ÿ Se 
ASL, Luni 
x soft shower 
me AP 10 


One of the consequences of the Gell-Mann (!) model of displaced charge 
multiplates for the strange particles is the prediction of a X-triplet (2*, 2°, 27) 
of which the charged counterparts (2*, 2") have been observed in various 
laboratories. The possible existence of the 2° has been inferred from the dyna- 
mics of the associated production of A°-0° pairs by FOWLER et al. (?), WALKER (3) 
in the diffusion cloud chamber and later by TRILLING and LEIGHTON (*) in 
the propane bubble chanber. Recently direct evidence for the existence of X° 


. GeELL-MANN: Suppl. Nuovo Cimento, 4, 848 (1956). 
. B. FowLER, R. P. SHUTT et al.: Phys. Rev., 91, 1237 (1953). 
. D. WALKER: Phys. Rev., 98, 1407 (1955). 
. H. TRILLING and R. B. LBIGHTON: Phys. Rev., 104, 1703 (1956). 
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and its decay into A° and y has been obtained by PLANO et al. (*) in a propane 
bubble chamber exposed to 1.15 GeV pion beam. In the three cases obtained 
by these authors the 2° was produced along with a 0°. 

We reproduce in Fig. 1 a cloud chamber picture taken at Darjeeling 
(2200 m) in which the evidence of the double production of a 2° along with 
a Kj, is as direct and convincing as those obtained by PLANO et al. (5) since 
we also observe the materialization of the photon arising from the decay of 
the X° over and above the decay of the resulting A° into a pion and a proton. 
In our case, however, the &° has been produced associated with a Ke whose 
decay products are also identified in the picture. 

A primary particle PI interacts at I in the lead plate (11.2 g/cm?) and 
produces the particle JO’ which stops and decays at 0’ into a muon O'S, 
(in the backward direction) of visible range greater than 92 g of Pb equivalent. 
The plates in the chamber are alternately 5.5 g of Cu and 11.2 g of Pb be- 
ginning and ending with Cu. Since the ionization of the particle O'S, in the 
last gap is appreciably greater than minimum we put its range R,= (97+ 5) g 
of Pb equivalent. This agrees remarkably well with the range of muons from 
K,, as obtained by the Ecole Polytechnic group and the G-stack collabor- 
ation (5). Moreover the range of the particle JO’ is (38 +3) g of Pb after 
it ionizes between 1.5 to 2 times minimum and therefore it is easily identified 
as a K-particle of mass between 900 to 1000 times the electron mass. Thus 
10'S, actually represents a K,, event in which the K-particle had a momentum 
of (369 + 10) MeV/c at the point of production. Taking the mass of Ki, ve 
be 966 m, its total energy is found to be (615 + 6) MeV at the time of emis- 
sion, where the statistical error shown is only due to the uncertainty in the 
total range. 

At O, just below the interaction, we find a V°-event whose secondaries OS, 
and OS, are definitely identified to be a pion and a proton from the following 
considerations. The secondary OS, has only a range of (15 + 6) g of Pb from 
the point of ejection and ionizes about three times minimum before stopping 
in 11.2 g of Pb. Hence it is surely an L-meson. Besides this, the particle 
OS, suffers a large angle scattering in the only plate (5.5 g of Cu) through 
which it passes and hence it is taken to be a pion. The other secondary OS, 
has a visible range greater than 66¢ of Pb equivalent and ionizes four to 
six times minimum in the last gap. Its range is more than 55 g of Pb after 
it ionizes more than two times minimum and this fact puts its mass much 
above the K-particle mass and very near to that of the proton. This secondary 
is therefore identified as a proton of total range (73 + 6)g of Pb at the time 


5 


(5) J. CRUSSARD, V. FOUCHÉ, G. Kayas, L. LEPRINCE-RINGUET, D: MORELLET, 
F. Renarp and J. TREMBLEY: Suppl. Nuovo Cimento, 4, 373 (1956); G-Srack CoL- 
LABORATION: Suppl. Nuovo Cimento, 4, 398 (1956). 
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of emission and has therefore a momentum of (709°) MeV/c and a total 
energy of (1174 +13) MeV. Having identified the decay products as a pion 
and a proton we take the V° event at O to be actually the decay of a A° whose 
momentum and energy are accurately determined from the total energies of 
its decay products and the known mass 1113 MeV of the A°. They are 
(774 + 28) MeV/e and (1356 +17) MeV respectively. As the secondaries OS, 
and OS; are observed in a plane almost in the line of sight of both the right 
and left views, the angle between them is difficult to measure. The difficulty 
is all the more enhanced by the large angle scattering of both the secondaries 
in the very first plate they 
pass through (see Fig. 2 C). 


soi vi 
Nevertheless the direction of 5, . 
the A° and its two secondaries P ner UN 
can be calculated from their i y sì 
known momenta and the dyn- } 

amics of a A® decay. We have a "e 

therefore determined as accura- 4 
tely as possible the point O, 

the direction of the pion OS, 2 
in space and the A°-decay plane , 5 
and then constructed the di- Fig. 2. - Schematic drawing of the whole sequence 
rection of the A° and proton of events of Fig. 1 with connected non-ionizing 


links (dashed lines), and various angles in space. 
« A » denotes the production of X° and K, along 
with the latter’s decay. « B» shows the disin- 
directions of the primary PI tegration of X° into A° and y with the latter’s 
and the K,,-particle IO’ do not materialization into four electrons a, b, €, d. 


change abruptly and their plane « C » depicts the decay of the A° and the large 
angle scatterings of both of its decay products. 


as required for a A°-decay. The 
result is shown in Fig. 20. The 


is fairly correctly determined. 
It is found in the stereo projec- 
tion that the apex 0 of the A°-decay does not lie in the plane determined by PI 
and JO’. In fact the true direction of the A° as determined previously makes an 
angle of (6-+2)° with the plane PIO’ pointing thereby to the possible emission of 
a > as an intermediate particle just as in the cases of the A°-0° pairs observed 
in the Brookhaven diffusion chamber and the Berkeley propane bubble chamber. 
An exactly similar double production of a A° and Kj, observed by BARKER 
et al. (5) does not however exhibit any non-coplanarity in the A°-K, plane. 
The range, momenta, and energies of the various secondaries and also the 
angles between them are given in Table I. 


(6) K. H. BARKER, M. S. Coates and B. R. FRENCH: Suppl. Nuovo Cimento, 4, 
319 (1956). 
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TABLE I. — Range, momenta and energies of the various secondaries and the various angles 
between them. 


Observed M ; Total tea 
Particle ech range in ei energy mee 2 Ù 
: tified as gJem? of Pb in MeV/c Mavi etwee 
IO' Ki 38+3 . 369 +10 615+ 6 WEST 
OG ut 9745 225+11 249-F10 | KE US AC ES 
OS, x 15+6 yee 182+11 | mp =63° +3 
È (calculated value) 
OS; p 73-46 Os 1174413 
leulated 2 
IO IN Selena 774428 1356+17 
JUS) Y 114+15 114415 | A0°-y =38° +4 
>a (calculated 867 +31 1470422 | nt-20 = 5.5°+1.5 
vela (calculated value) 
ILE Tt 1221 +38 1229 +40 mt-Ky= 13 +2 


We now come to the most interesting feature of the picture reproduced 
in Fig. 1. It is the presence of a photon initiated cascade of four electrons 
abed of which abe stop in 5.5g of Cu and «d» produces another pair of 
electrons stopping in the next plate (11.2 g of Pb). The apex of this cascade 
shower when joined in space to the point of interaction J gives a direction 
which is well within the core of the electron shower and represents, in our 
opinion, the direction of the photon emitted in the decay process, 


(1) 2° — AP +7... 


We shall now try to make a quantitative estimate of the energy of the 
ejected photon from (i) the dynamics of the decay mode (1) and also from (ii) 
the nature of the soft shower produced by it. Firstly we take the 4° and A° 
masses to be 1.187 GeV and 1.113 GeV respectively and conserving momenta 
and energies in (1) get the energy Æ, of the emitted photon as 


0.0847 


2 E, = 
©) Eno — pao Gos 0 


GeV... 


y 


where E, and py are the energy an momentum of the emitted A° in GeV 
and GeV/c and 0 the angle between the A° and the photon. The direction 
of the photon JS as determined by joining the apex of abcd to I in the 
stereo-projection, makes an angle of (384-4)° with the calculated direction 
of A° as depicted in Fig. 2B. This value of 0 combined with the known values 
of E, and p,. from Table I gives E,= (114 +15) MeV. Secondly the pho- 
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ton apparently traverses à distance of 2.3 radiation lengths (5.5 g of Cu and 
11.2 g of Pb) before giving a maximum number of four electrons and has there- 
fore and average energy of 125 MeV according to the cascade theory. Thus 
we find à remarkable agreement between the dynamically required value 
of the energy of the emitted photon and that obtained from the nature of 
the soft shower initiated by it according to the cascade theory. 

We therefore conclude that a X° was originally emitted at I along with 
the Kj and the X° immediately (in <10-® s) decayed into a A° and y. 
A schematie diagram of the whole event with connected non-ionizing links 
(as explained above) and various angles in space is given in Fig. 24, 2B and 20. 
Fig. 2A gives the production process together with the Ky decay event, Fig. 2B 
the decay of &° and the materialization of the emitted photon and Fig. 20 
the decay of A° and the large angle scatterings of the resulting pion and proton. 

Knowing the value of the energy of the emitted photon and also the energy 
of the A° from its decay products, the energy and consequently the momentum 
of the X° is obtained from its predicted mass (1187 MeV). So far we have 
not made any assumption as to the nature of the collision which has resulted 
in the associated pair and which has taken place in a lead nucleus. We notice 
that no other charged particle has emerged from the collision except the K,,. 
We therefore believe that a simple collision has occurred with a peripheral 
neutron and as a first approximation we balance the transverse momentum 
of the X° and K,, assuming the neutron to be stationary in the laboratory system. 
The angle which the K makes with the incident primary has been measured 
to be (13+2)° abd hence the angle of emission of the %° is found to be 
(5.5+1.5)°. Conserving longitudinal momentum, the momentum of the in- 
cident particle is obtained as (1221 +38) MeV/c. In order to conserve stran- 
geness and the number of baryons the incident particle has to be a positive 
pion so that the whole sequence of events is represented by the following 


equation: 


ie ee KS 


+ 
(3) A +y Dar 


soft shower 
ne = 10 


If we now balance the total energy on the two sides of equation (3), we find 
an unbalance of (83 +46) MeV which is only (44+2)% of the total energy of 
2085 MeV of the Z° and K} This unbalance is quite small and can be neg- 
lected. Even if it is real, it can be accounted for by the Fermi momentum 
of the interacting neutron inside the nucleus which may be as high as 
270 MeV/c for a lead nucleus. 


36 - Il Nuovo Cimento. 


1481 


554 M. S. SINHA and C. R. MALLIK 


We therefore conclude that this picture shows an event in which we find 
direct evidence for the production of an associated pair of 2° and K,, together 
with their subsequent decays as represented by equation (3). It is also in- 
teresting to note that if we make a relativistic transformation of the production 
angles to the centre of mass system, we find that the 2° has been emitted in: 
the forward direction and the K, in the backward direction. This agrees well 
with the observation of LEIGHTON and TRILLING (4) in the production © -K 
pairs. 


* Kk 
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RAVASS IS UIN On (+) 


Si è ottenuta una prova diretta della produzione associata di un X° con un Kee 
La caratteristica più interessante del quadro è che oltre a mostrare i decadimenti di 
un A°(+7x-+p) e di un Ki (> u*+v) mostra la materializzazione del fotone risultante 
dalla reazione 2° A++. I valori dell'energia del fotone emesso, quale risulta dalla 
dinamica del decadimento nonchè dalle energie degli elettroni molli prodotti si accor- 
dano assai bene. Si ricava il bilancio energetico totale dell’intero evento che è rappre- 
sentato dalla reazione seguente: 


wo DSK, 


1 N 
01 x o) 
Ÿ NX sciame molle 
ae 210 


(*) Traduzione a cura della Redazione. 
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D. B. FAIRLIE and J. C. POLKINGHORNE 


Tait Institute of Mathematical Physics - Edinburgh, Scotland 


(ricevuto il 17 Febbraio 1958) 


Summary. — The Low equation for a non-relativistic factorizable po- 
tential is investigated. The familiar solution for h(k), where h(k) is the 
scattermg amplitude, is obtained directly as a Low equation by pro- 
ceeding from the eigenstates of the total Hamiltonian to the eigenstates 
of the free Hamiltonian. The relation between the manifold of solutions 
of the Low equation for h(k) and the spectrum of the free Hamiltonian 
is immediate. The zeros of h(k) are the bound state energies of the free 
Hamiltonian. 


1. — Introduction. 


In a previous paper (1) we discussed the relation between the manifold 
of solutions of the Low equation, discovered by CASTILLEJO, DALITZ and 
Dyson (?), and the existence of «hidden structure », that is to say bound 
states of the free Hamiltonian. In particular the Lee-type model of Dyson (8) 
and NORTON and KLEIN (*) was discussed in detail. In this paper we wish to con- 
sider the simple case of a non-relativistic factorizable potential which has been 
previously discussed by HAAG (5). He showed how one can construct a Low 


(1) D. B. FarrLie and J. C. POLKINGHORNE: Nuovo Cimento, 8, 345 (1958). 
Referred to as I. 

(2) L. CastILLEJO, R. H. DaLITZ and F. J. Dyson: Phys. Rev., 101, 453 (1956). 

(3) F. J. Dyson: Phys. Rev., 106, 157 (1957). 

(4) R. E. Norton and A. KLEIN: Phys. Rev., 109, 584 (1958). 

(5) R. HaaG: Nuovo Cimento, 5, 203 (1957). 
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equation for this model starting from an expression for the eigenstates of the 
total Hamiltonian in terms of the eigenstates of the free Hamiltonian. The 
resulting equation for the scattering amplitude h(k) can be solved by the 
trick (5) of considering H(k) — 1/h(k) and obtaining from mathematical con- 
siderations the equation satisfied by H(k). The manifold of solutions arises 
from the fact that h(k) can be given an arbitrary number of zeros which pro- 
duce an arbitrary number of poles of H(k) and a consequent arbitrariness in 
the expression for H(k). 

In this paper we adopt a converse procedure and use an expression for the 
eigenstates of the free Hamiltonian in terms of the eigenstates of the total 
Hamiltonian. This yields a Low-type equation which turns out to be just 
the equation for H(k) previously obtained by indirect mathematical arguments. 
The way in which we obtain this equation makes it clear that for this model 
the poles of H(k), or equivalently the zeros of h(k), correspond to the bound 
states of the free Hamiltonian. 


2. — The reciprocal Low equation. 


We consider a non-relativistic system with H, the free Hamiltonian; 


+ 
H — H,+V the total Hamiltonian; |k) the set of eigenstates of H with eigen- 
values £, and satisfying the condition that 


(1) Lim exp [iH,t] exp [— iHt]|k>) =|). 


where |k> are eigenstates of H, with eigenvalues H#,. We shall suppose that 
in addition to the |k>, H, possesses a set of bound-state eigenstates |x) of 
energies H,. The complete set of eigenstates of H,, both |k> and |a), will 
be denoted by |!» and we have 


(2) DID 


The usual way of setting up a Low equation for this system is to consider 
the expression 


(3) | =| Vik), 


il 
H — E, +ie 


(5) G. F. CHEW and F. E. Low: Phys. Rev., 101, 1570 (1958). 
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which gives the eigenstates of H in terms of the eigenstates of H, and leads 
to the equation 


cd LALSZOANA dea 
JE — Hy — ie 


(4) VIE =-— I ||. 


If the potential V is factorizable so that we may write 


(5) UV |b = Au*(l')u(1), 


then (4) becomes the Low equation 


RENE [w(h") |? |A(k")|2, 
(6) eta ti 
where 4(k') is defined by 
ys = ee 
(7) <k iy ik) al u(k') 2 
t(k') 
8 h({k’) = —— 


The summation on the right hand side of (5) represents a combined sum and 
integral over the energy spectrum of H. 

We reverse this procedure and consider H, = H—Y. Then for an eigen- 
state |k> of H, belonging to the continuum we have 


+ + 
(9) |k> =|k> + Vib, 


H, — HE, +ie 


as the analogue of (3). Since the states |k> and > coincide at t=co the 
choice of the sign of ie in (9) is the correct one. Equation (9) leads to 


+ + 
5h + 3 <K|V|bA|VIk) 
kV oe <A Vba i =, 
(10) De UALS Dons er 
This equation is equation (16) of I derived from a different view point, and 
we shall now use it to solve the Low equation. From equations (5) and (7) 
and the identity 


M VIE = > (DAV | , 
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we deduce that 


(12) xe Dem = 7h 

We have the identity 

(18) GEV [BY = SE EDAV YEE , 
and from (5), (12) and (13) it follows that 

(4) dr = AT 


Substitution from (7), (8) and (14) into (10) yields 


1 1 1) |? 


The summation on the right hand side of (15) is over the energy spectrum 
of H,. It is an integral over the continuum states |k> and a sum over the 
bound states |x). That is 


à eee Le | (cc) |? 5 LIURIES 
CO nevi) tirata 


This is just the form of the solution of the Low equation (6) corresponding 
to the generalized R-function 


clue 
E,—E,— ie 


(17) R(k) = > 


in which H(k) has poles at H,= H.. 


3. — Discussion. 


For this model we have given a quasi-physical explanation of the trick 
used for solving the Low equation. This explanation has made it clear how 
the properties of the solution correspond to the spectral properties of H,, the 
free Hamiltonian. In fact the zeros of the scattering amplitude occur at the 
energies of the bound states of H,. 

There is always an intimate connection between the zeros of h(k) and the 
bound states of H,, though this relation is not always quite as direct as that 
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discussed here. In I we considered the Dyson model, which differs from the 
model treated here in that equation (5) is only true if one state is a continuum 
state and the other a bound state. There the relation between the « hidden 


structure » and the generalized R-function was of a different, but equally 
simple, type. 


ERAS Si os UN EEO mre) 


Si studia l’equazione di Low per un potenziale fattorizzabile non relativistico. La 
soluzione familiare per h(k), dove h(k) è l'ampiezza di scattering si ottiene direttamente 
come equazione di Low passando dagli autostati dell’hamiltoniana totale agli autostati 
dell’hamiltoniana libera. La relazione tra la famiglia di soluzioni dell’equazione di 
Low per h(k) e lo spettro dell’hamiltoniana libera è immediata. Gli zeri di h(k) sono 
le energie degli stati legati dell’hamiltoniana libera. 


(*) Traduzione a cura della Redazione. 
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(ricevuto il 21 Febbraio 1958) 


Summary. — The asymptotic expansions of the scattering amplitude 
at high energy obtained previously for a Schrodinger equation of motion, 
are extended to the cases of a Dirac and Klein-Gordon particle. The 
method followed here is based on the extension of the so called Gel’fand 
and Levitan operator which gives a uniform connection of two complete 
bases in the Hilbert space. A compact expression for the general term 
of the expansion is found. This allows its evaluation by means of re- 
currence relations. 


1. — Introduction. 


We wish to consider in this paper the same type of problem which we have 
already treated for the scattering of a non-relativistic scalar particle by an 
external field (1). The same restriction of a statical spherical symmetrical 
potential will be maintained, but we consider particles of spin 0 and 4 moving 
with velocity near to the limit value of the light velocity c. 

We deduce some expansions of the phase shifts relative to a fixed angular 
momentum and parity, in powers of (1 — 62), where 8 = v/c is the ratio of the 
particle’s velocity v and c. 

The mathematical tool used in order to arrive at the expansions is closely 
related to that of the Schrédinger case (1). Indeed, we start from an integral 
operator O which maps the waves ©, distorted by an external field V, to 


(1) M. VERDE: Nuovo Cimento, 6, 340 (1957). 
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the corresponding waves @ distorted by a different external field V. © has 
the peculiar feature of beeing independent of the energy. This is particularly 
suited to the study of the behaviour of ® for large E, when the corresponding 
behaviour of ®, is known. 

The same study done with the mapping of ®, into © by means of the 
familiar operator, which contains the appropriate Green’s function G of the 
Ds space, is much less powerful because G is energy dependent. 

The operator O corresponding to the cases treated in this paper has not 
been given before (?) and it is presented here without proof. Since its deduction 
requires some extended considerations, and it is especially important in con- 
nection with the inversion problem of providing the potential when the spectral 
function is assigned, it forms the subject of a separated investigation (*). 


2. — Potential scattering of a Dirac particle in a central field. 
The Dirac equation of motion for a partial wave of angular momentum 7 


and fixed parity reduces, as it is well known, to the following two component 
wave-equation 


220 lai: a 
(Ed) 10% o + Mos + Le o, + V (a) Lg (E, x) = Eg(E, x), 


where o, are the familiar Pauli matrices. 
It is convenient in the following to work with an equivalent form of Eq. (1.1) 
performing a unitary transformation by means of the operator 


de _ 


ta Da, ‘i 2 I 
U = exp [i0,0] with d V (a) 
In the new frame of reference Eq. (1.1) becomes: 
ied 

(1.2) — 102 = + 2) DE, x) = EDE, x), 
where 

er j+4 : 
(1.3) Q = exp [i0,0] (mos + 3 0,| exp [— ia,Q]. 


We consider now the motion in a different potential V,, which in part- 


(2) There is a tentative done for the Klein-Gordon equation of motion, which is 
however incomplete in our opinion. See E. CoRINALDESI: Nuovo Cimento, 11, 468 (1954). 
(3) M. VERDE: in course of publication. 
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icular may vanish. The Dirac equation corresponding to (1.2) will be 
Me 20 
(1.2') mr. + Q0| DO(E, x) = EDO(E, x). 


Choosing the ® and ®% so that they have the same behaviour for x= 0, 
namely 


544 


(1.4) D, DO ~ (rs ) + O(ai+8), 


the completness relations for ® and ®" will be written: 


famose, æ)D;(E, x") = Oxp(@ oe x") TUNIS x) D> (EL, a") ’ 


where 0(£) and o,(#) are the so called spectral functions. do(£), do,(E) have 
singularities of the type of the Dirac functions 6(# — E,), — m<E,<+ m, for 
the eigenvalues E, corresponding to ‘the bound states of Eq. (1.2) respect- 
ively Eq. (1.2). 

In complete analogy to the Schrédinger case, there is a connection between 
the two bases ® and D® of the Hilbert space, of the type 


(1.5) D(H, x) = (1 + iSo:) DO(E, 2), 


the matrix elements of S being 


+o 


(1.6) Sag(E) = fate 00) 


ce] 


D,(E'x) OD (E'2) 
TIP Da 


This follows from the Gel’fand and Levitan equation (*) equivalent to (1.2) t.i. 


+ © a 


(1.7) O(E, x) = DO(E, x) — Jane, n fo, y)BO(E", y) dy , 


ES 0 
and the continuity equation for D® namely 


— DE", 2) DE, x) = (E — Bf own y) P(E, y) dy . 


0 
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The form Eq. (1.6) of the operator S is well suited to the study of the behaviour 
of ® when HL is large, since we have to evaluate the moments 


K(x, 2) = |e d(o — 09) ®,(H, x) G9 (E, x), 


in those cases which allow an expansion in inverse powers of Æ of the same 
operator S. 

It is now possible to establish simple recurrence relations to which the 
moments have to satisfy and that serve also to their evaluation in terms of 
the difference of the potentials V(æ) — V,(x) =6V. Putting: 


+ o 


Ks6®, y) LT = 0) P(E, x)DE E, y) ’ 


= © 


and operating on it with io,(0/0x), from Eq. (1.2) one obtains 


CK (x, 
re ig, OY _ OK, y) — Ko(a, y) 


If one operates on the left of the same K”(x, y) with io,(0/0y), remembering 
Eq. (1.2’) we get also: 


OK (x, y) 


(1.81) ay 


16, Kt, yy QO(y) — Ka, y) . 


Since o, anticommute with 2 and 9% one deduces: 


AK (2,0) (aK aK 
de pores dx 


| = i(Q(a)o,K(w, a) — K(a, 2)o,Q(w)) + 
+ i(0,K%+D(x, à) — Ka, w)o,) . 


This is the desired recurrence relation, that can be put in an equivalent but 
more convenient form introducing the operators 


TE) = K(a, 2) +0, K(a, co, . 
One obtains then: 
(1.9) 


Ce ee 


Dio, 8° = ——— — (Qo, Pt’ — 0, TS 20). 
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Thus the knowledge of the first moment T° and the boundary condition 
for x= 0 allows the determination of all other higher moments. It is easy 
to deduce from the Eq. (1.7) where only Æ(x, y) appears, by a simple in- 
tegration by parts that 


(1.10) 10, 12 Oa) O(n)? 
The boundary conditions Eq. (1.4) give 
T'P(0) =— T(0)o; , 


which is sufficient for the integration of (1.9). One deduces then from (1.10) 
and (1.9) 


TOY (x) = cost = — T(0)o, = (27 +1)t0,(8V) pao 5 
EM ee = (2 — 2) — (j + 4)(8V).(2 — 2)io, , 
TO) = me [dV a Lj 4 v(| le de i | 


7 een a rare Oe 
av 


and so on. 

We have assumed here that (3V)' vanishes for x = 0, otherwise the moment 
T’? would diverge. 

For the evaluation of the scattering amplitude we need the values of the 
To = T™(co) when «=co. The first moments calculated above are: 


T® = m(exp [2i0,4] — exp [2i0,4])o, , 
TO = (25 +1)(8V oto, , 
TO = (j ARIANO te 


i 4 : SV) s 
PO = ne sA = (4 pf 2 dx + mo, exp [io, 5A] sin SA + 


+ (4G — HBV). + VO) TE, 


where 


A = [rinar, À0 [rar Sy ee Vhs KO). 
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Remembering the unitary transformation which brings the g’s into the Ds. 
we write 


yp = (exp [— to, 8Q] + i exp [— i0,Q]8 exp [ic,.Q]o,)y . 


Since at large distances 


~ 
| 
| 
| pole 


sin (er = x) 
1—o, E+ml+0; \ 2 
pe ~ = 


2 k 2 f i + 


j , 
cos | kr — 5 a) 


bole 


it is easy to derive for the tangent of the phase shift, the expansion 


k ( 1 Ê dA 
\ o Ci to i 3 A) 
(1.10) tg 0; E tg À i ! E2 in D 4 sin 2 a) | 
| 1 f(3vy 
DI 7 1\2 un 
tn sn) pie 


+0(%)}. 


The phase shift 0; corresponding to the same j but opposite parity in the 
same order of approximation results to be the same. 


3. — Potential scattering of particles with spin zero. 


If g(E, x) and g®(E, x) are the stationary eigenstates of a Klein-Gordon 
particle in a central potential V(x) respectively V,, so that 


INF) “oi 
{we vx mu x + aah? a 
(251) P 
+1 o? 
(Eva pu? — per = 0, 


and we choose the normalization g, yp ~ «+! when «= 0, there is again a 
uniform connection between the two complete bases of the Hilbert spaces 
expressed by the equation 

+ o 
cos (Q — 0) —] do(L’) 


AS 


(2.2)  p(E, a) =pE, 2) 


o(£", ap (E, a), 
RISE: i 


x (È n P(E', æ)p®(E",x) 
+ "(Bs 2) | dace pane 


— © 
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where the dashes indicate derivatives with respect to the radial variable x, 
and the other symbols are defined analogously to the spin } case. 
Namely 


Qe = [ave an do = de — de; 
and 
Paes + 0 
| do(E)p.(E, &)p,(E, #')(01)yp = Sap(@ — 2) =] do Eps (E, 2g (2, 1 (01),8 » 


(0) 


where p,, g; are two-components wave-functions 


p(E, x) PE, 0) 


4 
= ME) er PE, x) 


Eq. (2.2) is an immediate consequence of the Gel’fand-Levitan form of the 
Klein-Gordon equation. 
Indeed this form reads (8), 


(2.4) PAH, v) = = COS (= Q°)p PE , ©) — fut), 99 En dy 5) 


and the matrix elements of K are 


co 


(2.5) Kxp(, y) = rfactor , pp (E, y) . 


Eq. (2.2) follows from the continuity equation 


x 


(2.6) WB E) [9 Morn? XE, y) dy = 


0 


= @U(E, x) 


0po(E, x) = 0g (EH, x) 
Ca É 


after substituting Eqs. (2.5) and (2.6) in Eq. (2.4). 
We remark here that the form of the transformation given by Eq. (2.2) 
is very similar to that already established (1) in the case of non-relativistic 
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kinematics (cfr. Eq. (1.9) of the paper I), we have here a term cos (Q — Q®) 
instead of 1, as it must appear since the phase shifts approach a finite limit 
as Eco in the relativistic case. The moments 


+0 


K,,(@) =| do(£)E"o(E, x)p(H, x), 


— © 


Ha) = | do(E)E"p(2, aq (H, «) 


— © 


needed for the asymptotic expansions, are of course different from the non 

relativistic ones and in particular they satisfy to slightly more complicated 

recurrence relations due to the new form of the equation of motion Eq. (2.1). 
It is now straightforward to apply the same procedure followed in the 

non-relativistic case, although one is forced to more lengthy calculations. 
One arrives at the result 


reel 1 df oA \ Ub +41) SV) da 
2.7) te 6, =—— te | alex x 
ATO FA Ca z| si (; sin 2 gi "sin 254) @ dal e (Fa) 


If we compare this expansion with the corresponding one for a Dirac part- 
icle as given by Eq. (1.10) we find that they differ by a term 


k 1 (avy 
== == —_ = 4 a er 
pte maga (UD) —G+ nf n 


0 


which is entirely due to the spin-orbit coupling. Indeed {(1+1) —(j+4)? co- 
incides with the expectation values of the operator L-o in the substates 
j=1+4+ and j=1—34 and (1/x)(d/dæ)9V represents the familiar radial de- 
pendence of the Thomas precession. 

This result can be also obtained iterating the Dirac equation and retaining 
terms up to the order 1/E?. It is preferable to work directly with the Dirac 
equation, because of its linearity in the potential. This yields considerable 
simplicity in the actual calculations. 


4. — Final remarks. 
The expansions (1.10) and (2.7) serve mostly to follow quantitatively at 


large energies how the phase shifts change in correspondence to a different 
shape of the interacting potential. Naturally they are useless for large impact 


1495 


568 M. VERDE 


parameters. In this case, however, the incoming particle is moving outside 
the region where SV is appreciably different from zero and a perturbation 
expansion can instead serve the scope. It is probably also possible to obtain 
expansions valid for large angular momenta, following a procedure similar 
to that used in the case of non-relativistic kinematics (*). 

For the total scattering amplitude at high energies, since no difference 
exists between the partial waves corresponding to the same j and opposite 
parities, the same considerations already given in the Schrodinger case will 
continue to be valid (5). 


(4) See Sect. 3 of the paper mentioned in I. 
(5) See Sect. 4 of the same paper. 


RIASSUNTO 


Si ottengono sviluppi asintotici delle ampiezze d’urto per grandi valori dell’energia, 
in problemi di urto in un campo esterno a simmetria sferica, di particelle di Dirac e di 
Klein-Gordon. Viene stabilita una forma integrale per il generico coefficiente dello 
sviluppo. Essa ne permette la valutazione mediante formule ricorrenti. 
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à Pr = 


NDAZIONE FRANCESCO SOMAINI 
PRESSO IL TEMPIO VOLTIANO, A COMO 


Pa 


BANDI DI CONCORSI 
AL PREMIO E ALLA BORSA PER IL 1958 
Con lo scopo di premiare e incoraggiare nel nome di ALESSANDRO VotLra gli studi 


di Fisica in Italia, la « Fondazione Francesco Somaini », presso il Tempio Voltiano a 
Como, indice i seguenti concorsi. 


A) Concorso al “Premio Triennale per la Fisica Francesco Somaini” per il 1958. 
di L. 1500 000 (un milione e cinquecentomila) nette, da assegnarsi al concorrente che, 


fra quelli che la Commissione Giudicatrice giudicherà in senso assoluto meritevoli del | 
premio per i risultati conseguiti nello studio della Fisica durante il Triennio 1° Luglio 1955- __ 


30 Giugno 1958, sia, a parere della Commissione stessa, il più meritevole. 


B) Concorso alla Borsa Francesco Somaini per lo studio della Fisica ’’ per il 1958. 
di L. 750 000 (settecentocinquantamila) nette, da assegnarsi al concorrente che, tra 
quelli che la Commissione Giudicatrice giudicherà in senso assoluto meritevoli della 


Borsa, verrà dalla Commissione stessa giudicato il più meritevole, sia per titoli, pre- 


parazione scientifica, lavori già svolti e risultati già conseguiti nella Fisica, sia anche 
per il vantaggio che gli studi, per i quali è richiesta la Borsa, possono portare allo svi- 
luppo della Fisica, in Italia. 


1. — Ad entrambi i Concorsi possono prendere parte singolarmente i cittadini, d’ambo 


i sessi, italiani e svizzeri del Canton Ticino purchè di stirpe italiana. Sono esclusi dal — 


Concorso i membri della Commissione Amministratrice e della Commissione Scien- 
tifica della « Fondazione Francesco Somaini ». 


2. — Le norme particolareggiate dei singoli Concorsi verranno pubblicate in appo- 
sito volantino che potrà essere richiesto dagli interessati alla Segreteria della Fonda- 
zione presso il Tempio Voltiano a Como. 


8. — La domanda, i documenti, i lavori, ecc., presentati dai singoli concorrenti — 


dovranno pervenire, tra il 1° Gennaio e le ore 12 del 1° luglio 1958, alla Commissione Am- 
ministratrice della « Fondazione Francesco Somaini » a Como presso il Tempio Voltiano. 


4, — Il Premio Triennale per la Fisica potrà essere anche conferito a uno studioso 
che non abbia preso parte al Concorso, ma sia stato segnalato da un Membro della 
Commissione Giudicatrice, con proposta motivata, come meritevole di particolare con- 
siderazione oppure ritenuto degno di premio dalla Commissione Giudicatrice, indipen- 


dentemente da ogni segnalazione. 
* OK OK 


La procedura dei suddetti Concorsi è regolata secondo lo Statuto della Fondazione, 
il quale è ostensibile a Como presso il Tempio Voltiano ed è depositato negli atti del 
Notaio Dr. Raoul Luzzani di Como. 


Como, dal Tempio Voltiano, il giorno 6 Agosto 1955. 
Il Segretario Conservatore del Tempio Il Presidente Sindaco di Como 


CESARE MORLACCHI PAOLO PIADENI 


[Il bando in esteso di detti Concorsi è stato pubblicato nel fascicolo di Gennaio del Nuovo Cimento 1956], 
x, 


Gi FASCICULES SPÉCIAUX 


DU 


SUPPLEMENTO 


AUX VOLUMES DES SERIES 1X) ht x 


3 NUOVO CIMENTO 


qu. e Comptes rendus de Colloques, Congrès, Conférences internationaux, orga- 
KR nisés par la Société Italienne de Physique. 

fo” Colloque International de Mécanique Statistique (Florence, 1949). 

Es 1° Congrès International sur les Rayons Cosmiques (Côme, 1949). 


| DSG Colloque International sur les Ultrasons (Rome, 1950). 


Ch 

«Sr Colloque International d’Optique et Microondes (Milan, 1952). 

TR Congrès International sur les particules instables lourdes et sur les événements 
13 de haute énergie dans les Rayons Cosmiques (Padoue. 1954). 

a) Colloque International sur la Physique des solides et liquides (Varenna, 1954). 
LME Colloque International d’études sur VInfrarouge (Parme, 1954). 

wu; Conférence Internationale sur les particules élémentaires (Pise, 1955). 

AO Colloque International des études sur Vionosphére (Venise, 1955). 
i e: i __Colloque International sur les constantes fondamentales dela Physique (Turin,1956). 
“Wa Colloque International sur les rayons cosmiques (Varenna, 1957)*. 

oe: Colloque International sur l’état condensé de sistèmes simples (Varenna, 1957)*. 
"RI ART. Colloque International sur les mésons et les particules récemment découvertes 
Di è : (Padoue-Venise, 1957) (s’adresser à V Istituto di fisica dell’ Universita - 
Fo; «Padova - via Marzolo, 8). 

STI * sous presse 


_ —@ Cours tenus à Varenna à l’école internationale de Physique de la Société 

Italienne de Physique. 

Révélation des particuls élémentaires, avec référence spéciale à la Radiation 
cosmique (Varenna, 1953). 

È Physique des particules élémentaires et machines accéleratrices (Varenna, 1954). 

> Questions de Physique Nucléaire: structure et modèles du noyau, nouvelles espèces 

asi d’atomes (positronium, atomes mésiques), Optique neutronique, moments 
nucléaires, processus physiques du réacteur nucléaire (Varenna, 1955). 

| Propriétés magnétiques de la matiére (Varenna, 1956). 

CA Physique de l’état solide (Varenna, 1957)*. 


* sous presse 


@ Rapports sur des questions diverses de Physique. 
Premier, Deuxième, Troisième, Quatrième, Cinquième Rapport sur les travaux 
de Physique en Russie (1953, 1956). 
= MORPURGO et FRANZINETTI. « An Introduction to the Physics of the new part- 
STE icles » (1957). 


Pour les commandes s'adresser 5 
NICOLA ZANICHELLI EDITORE 
34, Via Irnerio - BOLOGNA 


ra 


IL NUOVO CIMENTO Vou. VIII, N. 4 16 Maggio 1958 


Photoproduction of 7--Mesons Near Threshold 
and the Value of the Panofsky Ratio (*). 


A. BALDIN 


Physical Institute of the USSR Academy of Sciences 
Department of Mathematical Physics, University of Birmingham - Birmingham, England 


(ricevuto il 24 Febbraio 1958) : 


Summary. — An examination of the cross-sections for the reactions 
2n + nt 
DI + d > 
5 | 2p + n 


near threshold was performed. The comparison of the results of this 
examination with experiment is presented, and it is shown that the 
negative to positive ratio is constant in the photon energy range 
(153 +193) MeV and equal to 1.4 (to an accuracy of 10%) in agreement 
with field theoretic results. The disagreement of these results with the 
values of the Panofsky ratio and the S-wave scattering phase shifts is 
discussed from the point of view of the two-component r°-meson field 
hypothesis. A direct experimental test of this hypothesis is proposed. 


1. — Introduction. 


The interaction of low-energy z-mesons with nucleons is one of the most 
interesting problems in pion physics. The understanding of the interaction 
in the S-state leaves much room for improvement. The experimental data 
on this interaction are discussed in numerous papers and conference reports. 


(*) A part of this work was performed in the Physical Institute of the USSR 
Academy of Sciences and reported in the Padua-Venice Conference (September 1957). 
The other part of the work was performed in the Department of Mathematical Physics, 
University of Birmingham where the author is at present a visitor. 


37 - Il Nuovo Cimento. 
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At the last conference in Rochester it was stated (*?) that the different exper- 
imental data on this interaction (the S-wave scattering phase shifts, the 
Panofsky ratio and the cross-sections for the photoproduction of 77 and x* 
mesons near threshold) are consistent between themselves within the limits 
of experimental error but inconsistent with the predictions from the threshold 
theorems and the dispersion relation approach (*). The threshold value of 
the negative to positive pion ratio in photoproduction o-/o* given by the field 
theoretical approach is equal to 1.4 if one uses for the renormalized coupling 
constant the commonly accepted value f? = 0.08 (*). From the experimental 
data on the Panofsky ratio, S-wave scattering phase shifts and o* (*) one can 
get for o-/o+ a value of about 2. 

Direct measurement of o—/o+ with extrapolation to the threshold (*) gives 1.8 
The discussion of this puzzle is the purpose of this paper. The direct measu- 
rements of the value of o-/o+ near threshold are not free from objection. These 
measurements involve the detection of the pions of a given energy in the 
reactions: 


(1) (2n+nt, 
Rai 
(2) [) 2p CETTE 


The measurement of low-energy pions is not a measurement of the threshold. 
cross-sections o- and ot because of the three body problem in the final state 
and the spread of the bremsstrahlung spectrum. One must have more detailed 
information on the reaction (2) before coming to conclusion about o— for the 
free nucleon. Then it is necessary to analyse in detail the mechanism of 
photoproduction of pions from deuterons taking into account all effects which 
contribute more than 10% (the accuracy of ot measurement). The theory 
which one can get as a result of such a discussion must be compared in detail 
with experiment because it is at present impossible to treat the problem without. 
making assumptions. 

In Sect. 2 we discuss the theory of reactions (1) and (2) near threshold. 
In Sect. 3 we compare this theory with the detailed experimental data (5) 


(1) J. M. CasseLs: Proc. VII Ann. Rochester Conference (1957). 

(?) G. F. CHew: Proc. VII Ann. Rochester Conference (1957). 

(3) Cfr. G. F. Coew, M. L. GOLDBERGER, F. E. Low and Y. Nampu: Phys. 
Rev., 106, 1345 (1957). 

(4) M. BENEVENTANO, G. BERNARDINI, D. CARLSON-LEE, G. STOPPINI and L. TAU: 
Nuovo Cimento, 4, 323 (1956). 

(*) Strictly speaking the ratio o/s does not depend on f?. However according 
to (*) it depends on unknown parameter N°. If one obtains N‘ from comparison 
ere vee experimental value of it, assuming f?= 0.08, the value of the ratio. 
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on reaction (2). As a result of this comparison we give the value and energy 
dependence of the squared modulus of the amplitude of the reaction Y+n > 
>p-+n- (Fig. 5). Our results are in agreement with the field theoretical 
approach and in disagreement with the experimental data on the S-wave 
phase shift and the Panofsky ratio. In Sect. 4 we briefly discuss the new 
situation in the pion-nucleon S-wave interaction (*). 


2. — Theory. 


The theory of the photoproduction of mesons from deuterons in the im- 
pulse approximation was discussed in several papers (°7). Our task is to 
investigate carefully the effect of the final state interaction in reactions (1) 
and (2). Some aspects of this problem we have discussed in reference (7). 

The total Hamiltonian of the system has the form: 


(3) H=H,+H'+V+V,, 


H, is the sum of the Hamiltonians for the free meson field and the interacting 
nucleons, H' contains the terms in the Hamiltonian involving the interaction 
between the nucleons. The meson field and electromagnetic field V and V, 
are nuclear and Coulomb potentials of the meson in the presence of the nucleons. 
The transition operator for the problem can be written in the form: 


il 
| | Ta = eA IEE I 
(4) TE" Mero mas a 


if one treats the last three terms of (3) as a perturbation. The first term of (4) 
is the usual impulse approximation. The operator H' has the form: 


(5) > exp life — k)t,] [(oK") + L’], 


v 


here the o” are the Pauli matrices, x and k are the momentum vectors of photon 
and meson respectively, K” and L” are the «spin flip» and «non-spin flip » 


(5) M. I. Apamovicu, V. I. VEKSLER, G. V. KUSMICHEVA, V. G. LARIONOVA and 
S. CHARLAMOV: Proc. CERN Symposium (Geneva, 1956), and private communication. 

The author is indebted to this sroup of experimentalists and especially to M. I. 
ADAMOvICH for detailed information about their results. 

(*) A more complete discussion will be given in a separate paper. 

(6) G. F. CHew and N. W. Lewis: Phys. Rev., 84, 779 (1951). 

(7) A. M. Batprn: Thesis (Physical Institute of the USSR Academy of Sciences, 
1953); CERN Symposium Proc., 2, 272 (1956). 
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parts of the transition operator for free nucleons. A simple calculation (’) 
results in the following expression for the cross-section of the reactions (1) 


and (2) (*) 


where p=(p,—p,)/2 is the relative momentum of the two nucleons and 
q = (Pi + p:)/2 


16 


Jk 
(270)? ie 


(6) lo 


2008 
Ial}:3 + gl "n cnr] SE — E;) dp dq , 


Il = fa exp [igr]pa dr and Ila = fo» exp [igr]pa dr , 


Pa = VEE AA deuteron wave function (*). 
| TT r 


The suffix «a» or «s » means that the co-ordinate part of the wave function 
for the relative motion of the nucleons is antisymmetrical or symmetrical 
respectively g,, and g,, were taken to be: 


ik al ; k 
(TT 3 Col (exp [ipr] — exp [—ipr]) , 


si FN. a : 1 — exp [— 276] exp |— iprl 
Cal exp [ipr] + exp [— ipr] 4 7 x 


(8) Pr = 


5) 


Boi 


for reaction (1). Here 6 is the phase shift for neutron-neutron S-wave scat- 
tering. This last wave function substantially exceeds the true wave function 
at small distance 7. The contribution from small r involves two other effects 
which are difficult to evaluate at present. The first is the competing process 
y+d — p+n and the second is nuclear multiple scattering of the meson. Mul- 
tiple scattering corrections give contributions of the order of — a/r, where a 
is the scattering length. For non-relativistic mesons a/r < 0.1 ifr = 1. Both 
of the above mentioned effects decrease the wave function at small distances. 

To take into account all these corrections let us assume that g,,, equals 
zero for r < 1. Such an assumption is also suggested by the A? dependence 


(*) Here we have neglected recoil terms in A and L depending on the velocity 
of the nucleons (7). They can give a contribution to our results of about TOS: Tie 
if one neglects such terms in discussing the m+ output from a deuterium target an 
error of about 20% can be involved. (This term can give contributions of opposite 
sign to the m+ and x amplitudes. The terms from the second part of (4) have also 
been omitted. For a discussion of these terms see below). 

(*) À = ¢ = w= 1 where u is the meson mass. 


dt 
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of photoproduction of mesons from nuclei. We therefore subtract from J,, 
the values 


1 
Olas = fr ar as EXP [igr ]padQ, . 
0 


To be on the safe side we should consider only those values of our variables p 
and q for which this correction to the cross-section is less than 10%. However, 
as we shall see later, this assumption gives rather good agreement with a pro- 
nounced experimental feature even in those regions of p and q where this 
correction amounts to about 25%. 

Evaluation of the integrals involved results in the following expressions 
for the cross-section of reaction (1) 


do 


6 2 + Bp, QL 
(9) SPE Gi g)|E|}+B(p,g)|L}}, 
where 
Boe p? + @\ | 4pq 1 1 x + (p +4) 
2 = p|1l = i 
ANP) a) A Mx Wee red apg Sri rap 
2 sin? 6 2Qqa pq ete o) s(n? Oe PG) al | 
= arete z - 1.8 — | In —— e 
3 | “erp IT 4\ a + (p—a)? 
1 sin 26/ 2qa pa ctg dì), e +(p+@?] 
a or (avete api 1.84 21 In ea (ge a dif , 
i 3a DISE a) 4pq 
NCIS oat l(a? + p? + 42)? — 4p2q? 
1 a2 + (p + 91 


I 
et pte a+ (p_-9)) 


A and B are practically independent of x. The region of the variables con- 
sistent with the conservation laws is 


Il 1 
0 <P < Va “ dx € (dq —)?] ; 


here q,=%x/(2-+(x/M)), M is the neutron mass e includes the binding energy 
of the deuteron and the neutron proton mass difference. A(p,q) and B(p, q) 
are peaked at p=—q. This peaking is connected with the weak binding of 
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the deuteron. For free nucleons the dependence of the transition amplitude 
on p and g would contain a factor $(p + 4q). 

In the region p< 0.2 A(p, 9) > B(p, q) because A(p, g)|K| is the spin 
flipping part of the cross-section and involves enhancement due to the final 
state interaction in the !S,-state of the neutrons. 

The correction 8Z,, gives contributions < 10% in the region p > q and 
about 25% in the region of small p. These corrections affect the nucleon 
S-wave contribution and have very little effect on other partial waves which 
give the main contribution at p > q and a small contribution for small p. 

For the discussion of reaction (2) it is necessary to take into account the 
Coulomb forces between all three particles. Roughly speaking the contri- 
bution of the Coulomb correction is of the order of magnitude of — z(e?/V) 
(V, relative velocity of interacting particles). The protons, being heavy part- 
icles, have much smaller velocities than mesons and must be more strongly 
influenced by the Coulomb interaction. But at the same time the nucleons 
have shorter wave lenghts and for them the estimate z(e?/V) is applicable 
only in the region of very small p. It is necessary also to keep in mind that 
the region of the small relative velocities (small p) is the region of enhancement 
due to the nuclear final state interaction. It is clear from this brief discussion 
that for the evaluation of the role of the Coulomb repulsion detailed calcul- 
ations are necessary. 

Instead of the functions (7) and (8) we take the exact Coulomb wave 
function: 


"I ù 1 TT 
(10) Ca) exp [ipr] > Co) exp 5 nl 

-|1( + in)|exp [ipr] Pin, 1, — i(pr + pr)], 
n= (e/V), F, hypergeometrical function 


1 — exp [— 210] exp [— ipr] 
ip P 


(11) 


> 2 exp [— 10,](1 — exp [— 2%6,])- 


exp [— ipr]W;(— in + 1, 2, 2ipr), 


TT 
-exp|—3y 


where W,, Whittaker function, dò, = arg 1(1+in); 6,—S phase-shift in p-p 
scattering. 


The method of evaluation of the integrals: 


(12) poe fer [inr] 22 = DR, na 
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can be found in reference (*). Expressing F as a contour integral (see Fig. 1) 
1 eat 

13 F(— in, 1, B) = — exp [i == ana 

(13) (— in. 1, B) gi EP [ip | P (1 i) ‘exp [ix] 


and changing the order of the ordinary and the 
contour integration one can get after calculation 


Arla? — p? + q? — 2iap]" 


(14 dl rao 
ee eS prose a 


The integral with Whittaker function cannot be 
evaluated in terms of elementary functions. Express- 
ing W as a contour integral (see Fig. 2) 


À a | 1 1\-*" exp [xz] de 
(15) Wi(în + 1,2, 2) = dari cp (1 ei i) & a? 


vi 


one can get after performing simple transformations 


x 1_-exp[2an] 


(16) [exp [ipr]W,(iy + 1, 2, — 2ipr) exp [inr] = 
Al) 


pq 270] 
L 
1 i 
-fexp [iy Int . - : 
Jes tin ETES ptqtia,, E GUT, 
i ptq—ia | peligna i 
27m 
Ca 
© exp [277] — sg 
£ 
i f, and f, are real functions which can be evaluated 
numerically. For small 7 one can use the expansion 
rata Sg INES sae 


© and f; can be evaluated with rather low accuracy. Expanding J, and 


taking into account only the first term in 7 permits simple calculation of 


integrals: 
| lasts MORO, 


(8) A. SOMMERFELD: Atombau, II (Braunschweig, 1956). 
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The results of these calculations are 


(17) se AP, 4) EA) I 


Mx 


4pq 1 1 net et a 
(08 prt gt) —4ptg? NS xp Pgo (pq) 


9 sj j 9 sj i 9 
pe e e, =(1 by arctg = — ) 
3 O*pq 3 pq CHENE TT 


x e 


3 2 2 2ap 
(18a) <A*(p,q) = = p (1 pra | C? I + 2n arctg = oe) 


[fz cos (6, + do) + fi sin (6, + 6,)]In 


ao L 2 _—_ m2 
(185) B(p, 9) = (i — 2n arctg es Bp, 4) ; 


and 


Ar 


DIL 


In the region of p < 0.1, A° < 3A and the Coulomb corrections are very 
important. The large difference between B° and B is not essential because B° 
and B are very small in this region. In the region of 0.2 < DIS 0 ATA 
and B° are smaller respectively than A and B by about 10%. At p = 0.4 the 
difference between 4° and A and B° and B is beyond the accuracy of our 
calculations. The numerical values for A, A° are given in tables I and II for 
the most essential region of our variables. 


TABLE I. - A(p, q). 


| | 
q 0.2 0.4 | 0.6 | 0.8 | 
P La | | | 
| | 
0.2 0.18 0.25 | 0.20 | 0.14 
0.4 0.13 0.30 | 0.26 | 0.16 
0.6 0.20 | 0.38 | 0.27 
0.8 0.09 | 0.30 0.49 
| | 
1 0.08 | 0.20 | 


(*) We do not expand (? here because one can show a eqs. (18) are applicable 
in the region of small p. At p < 0.1, f, & — x and fhe 
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TABLE II (*). — A°(p, q). 


q 0.2 0.4 0.6 0.8 
Pp 
One 0.16 0.22 ONE 0.12 
0.4 0.12 0.28 0.23 | 0.14 
0.6 0.20 0.38 027 
0.8 0.09 0.30 0.49 
1 | 0.05 0.15 | 0.37 
(*) A° is given without the correction from equation (24) because this correction is dependent 
on x and it is more convenient to take into account the factor (24) separately. 


To take into account the nucleon-meson Coulomb correction (*) it is 
necessary, according to (4), to multiply the transition matrix <p’, q', kjm,|T |i 
by the matrix 


(19) <P» q; k, 1 m (Vy Vs) 


Il 
Hey lis uN 
ii Pd k'j'm;», 


wheres Vi a/r Mie hr ret 

Ti, Do, T3, are the radius vectors of the protons and the meson respect- 
ively, 7 and m, the spin of the nucleon system and its projection. Straight- 
forward calculations give for (19) (+) 


Are? +2 


Cake + Be] 
: k— k’ 
Jesse COS (n ou Jens re 


Denoting the above integral by I pp (E — k') and introducing new variables, 
one can obtain the following equation for the correction of the transition matrix: 


ef, (yt dp’ dt al 
(21) 2 [rl 9 | en ales AP bs ga pa 
pa 2 Jp°—p? y e ER 
Î 


Me E 


(20) 6,0 


jj 0 


2q +2q'— k+ k') 


(*) This question was also investigated by M. MoravosiK. The author is indebted 
to Dr. Moravesik for informing him about this work and sending preprints of other 
works on photoproduction. 

(+) Here Bis the real physical cut-off taking into consideration inner bremsstrahlung. 
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Equation (21) can give a contribution greater than 10% only for meson mo- 
mentum less than — 0.3. A small momentum of the meson corresponds to 
large p (according to the conservation laws) and, as was shown above, small 
contribution of the nucleon S-wave part. So one can neglect the S-wave part 
of y,, and ,,, and I,, reduces to a combination of the 6-functions. 
Performing the integrations on the assumption that the integrand has no 
singularities except those written explicitly, one can obtain 
age ee peu) | 
p | 


(22) aoe ss (i | 
Uk —(q + p)/M] 

1 |k—(q— p)jM ; 
a ia air. 


For k < 0.3 it is possible to neglect the interaction between the nucleons. 
This means that <p, q, k|7|i> is a real quantity and it is necessary to take 
into account only the real part of the correction. For f +0 we have: 


Te? | il 1 


2 \ ik = (q- 


|k — (q + p)/M| | L= =) <P) 4 k|T|è}. 


pIM \ 
(23 ie 
Sai (lana) <2 
almost in the whole region of interest. 
Neglecting p?/(M?|k — q/M |?) we can obtain that the correction on meson- 
proton Coulomb interaction to the cross-section of reaction (2) reduces to 
the factor: 


2re? | 27e? 
24 (1 + — — | ~ (1 = = 
124) lk —q/MI., [k—x/2M | 


(1 | 0.046 | 
zi I = = == —— = — = — = = = . 
V (x2 — 1.09) — 0.32(x — 1)q? — 0.32xp? 


Our final result is practically equivalent to the operation of multiplying the 
cross-section by |¥(0)|?, where Y is the wave function of the meson in the 
Coulomb field of the protons. This means that it is valid only when the wave 
lengths of meson-nucleons relative motions are large in comparison with the 
Size of the system: 
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oO 
+I 
© 


For À — 1/x the correction is about 13%. For 1/À > « the correction (24) is 
not valid, but the real correction is beyond the accuracy of our calculations 
and must be put equal to zero. The last statement follows from a simple ana- 
lysis of equation (18b) which describes a very similar problem. 

The restriction (23) for the applicability of (24) practically coincides with 
the criterion of applicability of perturbation theory because in the region of 
interest p <1. 


3. — Comparison with experiment (*). 


For the most part our discussion deals with experimental results obtained 
in the Physical Institute of the U.S.S.R. Academy of Sciences (°). This exper- 
iment gives complete information about the process (2) for energies of the 
photons from the threshold up to ~ 190 MeV because all three final particles 
are detected. The coefficient 6 in the region 7 < 0.3 is negligible. This per- 
mits us to write 


fu a nay! dp Hi (p, q) dp 
i do Bary 
9 3 dpd 
TE RIT LE [fra ;(p, q) dp dq 


Comparison of this value with experiment gives |L|?/|K} < 0.1. This sup- 
ports the well known result that the S-state of meson and nucleon gives the 
main contribution to the cross-section of photoproduction of pions near thre- 
shold. Within the limits ~ 0. 

To check our theory the following characteristics of the process were ob- 
tained 


° do p? 
— = Ii T met 
1) la ag dg.) where nà 
and 
2) (a dp = g;(€) where Ee — (p —Q). 


These characteristics are very sensitive to the description of the relative motion 


(*) The detailed discussion of the experimental data will be published in ‘ie Journ. 
Exp. Theor. Phys. USSR. 


1507 


580 A. BALDIN 


of the protons. Experimental data for g,(7) and ¢,(e) are shown in Fig. 3 
and Fig. 4. The solid curves are the results of numerical integration in (26 
and (27) with do/dp dq from (17), (18) and (23) and with K |? =0.8-10-*" em?. 
The last value corresponds to the field theoretical prediction with the 
coupling constant f? = 0.08. 
As one can see the agreement with experiment is very good. This permits 
us to discuss the dependence of |K | on x. 


Ag (7) 


siena is 
002 004 006 008 000 O12 O1 006 


Fig. 3. Fig. 4. 


Experimental data for the total cross-section of reaction (2) are: 


| | 
x (in MeV) | 156.5 163.5 | 170 181 | 
Gal0 eme: 4.6 + 0.6 gh 2b OLR 9.0 + 1.0 | ik Le sy | 


Dividing o by the results of the numerical integration | [A q) dp dq 
one can obtain the dependence of | K |? on o presented in Fig. 5. In Fig. 5 
are also given data about |A“?|? from the im- 
portant work of BENEVENTANO et al. (4). It is 
necessary to stress that our results are not in 


n 
1K? 107? cm?) 


10 

og Nes eee disagreement with experimental data on o7/ot 
aS Fa n Gul of reference (‘), but only with the interpreta- 
be tion of these data. In the work (‘) the measure- 


dn ica x (MeV) È Ta. ST Pi i Te i @ 
ey eo 70 7 ment: of the o-/ot ratio was performed by 


counting only = and x* mesons of a given 
energy. As was shown above, this ratio is de- 
pendent on the energy of the meson. Coulomb 
interactions are responsible for increasing o- with decreasing meson energy. 
For high meson energy the positive Coulomb correction (meson-proton inter- 
action) is negligible but the negative Coulomb correction (p-p interaction) 
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can give an appreciable contribution. For low meson energy (< 10 MeV) the 
negative Coulomb correction is negligible but the positive correction can give 
a contribution of more than 10%. In the work of BENEVENTANO et al (4) 
such behaviour of o7/o* was extrapolated to very low meson energy and it 
was assumed that the resulting value 1.87 is the threshold ratio 6-/o+ for free 
nucleons (*). This result is in disagreement with our’s, but the highest value 
of c_/5* which was actually observed in the experiment of BENEVENTANO et al. 
was 1.5 at a meson energy of about (10--15) MeV. For such meson energies 
the positive Coulomb correction is, according to (24), about 10% and the 
negative correction, according to Table II and the energy conservation law, 
is negligible. Introducing this correction, we obtain excellent agreement with 
our result. Similar results can be obtained from the values of o~/o+ found 
in reference (*) and attributed to another «effective y-ray energy », but, as 
was pointed out in Sect. 1, such an analysis is not very conclusive. For the 
function 9,(7) and @,(e) the contribution of the positive Coulomb correction 
is very small. 


4. — Discussion. 


The situation can be summarized as follows: The analysis of experiments 
involves careful measurement of reaction (2) very near the threshold and shows 
that the o-/o* ratio for free nucleons is independent of the energy of the 
photon and equals 1.4 (~10% accuracy). This behaviour and the value of 
the o-/ot ratio is in very close agreement with the predictions of field theory 
with a coupling constant of f? = 0.08. 

These conclusions contradict only the interpretation of the previous exper- 
imental data on direct measurements of the o-/o+ ratio but not the experimental 
results themselves. The disagreement of our results with the data on the 
Panofsky ratio and S-wave phase shifts in pion-nucleon scattering is, some- 
what of a puzzle. From the o obtained in the present work and from 
the well known Orear data on |ò, — ò;3| one can obtain a Panofsky ratio of 
2.4+0.4 against 1.5 +0.1 as reported by CASSELS (1). If this disagree- 
ment is not a result of a trivial overestimation of the experimental accuracy, 
then it is necessary to look for a more speculative explanation. 

Let us suppose that there exists (apart from the commonly accepted 
r°-meson) a x$-meson which has isotopic spin zero. The first consequence of 
the existence of such a meson is that it is necessary to perform a new analysis 


(*) Some objection against the procedure of extrapolation used by BENEVENTANO 
et al. was also discussed by M. MoravosIK (preprint). 
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of experiments on r--proton interactions. It is then clear that the value of 6, 
can be changed. 

The second consequence is that if one assumes the mass of ordinary 7° to 
be approximately equal to the mass of the x -mesons, then the Panofsky ratio 
involves mainly the 7° meson and the observed mass 264 is the mass of 7). 
This assumption may be made plausible on the basis of the isotopic invariance: 
The values of (m*— m°)/m® (where m* and m° are the masses of the charged 
and neutral particles involved in the isotopic multiplets) are: for nucleons 
0.15% and for X-particles 0.3%. But for m-mesons it is 3.3%, i.e. about 
10 times larger. The last value is closer to the ratio (M, — M,)/My= 6.5% 
(where M, and M, are the masses of the X and A particles) than to the first two. 

The 7°-meson properties obtained from absorption of mesons in hydrogen 
and deuterium are now, according to our hypothesis, the properties of the 
Ty-meson, i.e. 7 is pseudoscalar and strongly coupled with nucleons. 

Then we see that the existence of x, would completely change the very 
much discussed connection between important low energy parameters in pion 
physics. I do not know of any experimental result which would be in direct 
contradiction with the assumption about the existence of 7}. Our knowledge 
of the transformation properties of the z-mesic field is based on numerous 
tests of the consequences of the isotopic invariance hypothesis. The only 
direct test of the transformation properties of the r° field (investigation of 
the reaction Y+d ->d+-x°) was discussed in earlier works (°°), but as was 
mentioned in (°) this test could not exclude the possibility of descritpion of n° 
by two wave functions y, and q,) (*). 

I should like to point out a direct experimental test of the existence of the 
second component 7° field. The reaction 


d + d — ‘He + 70 


is trivially forbidden for ordinary mesons by conservation of isotopic spin and 


allowed for the x. 


The existence of 7) can give the solution to one of the most interesting 
problems of the physics of elementary particles, the problem of the S-wave 


(9) A. BaLDIN and V. MikcHafLov: Dokl. Akad. Nauk SSSR, 91, 479 (1953). 

(10) A. BALDIN: Suppl. Nuovo Cimento, 3, 4 (1956). 

(*) Professor R. E. Prrpris kindly communicated to the author that Y. YAma- 
GUCHI has also pointed out that at present there exists no experimental data incon- 
sistent with the hypothesis of an isotopic scalar meson. (After writing this paper I 
have been shown a manuscript of Y. YAMAGUCHI in which he discusses some conse- 
quences of different assumptions about the nature of the x9-meson, but there is little 
overlap between our hypotheses). 
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pion-nucleon interaction. A discussion of this and other consequences of the 
existence of the 7° will be given in the next paper. 


It is a pleasure to acknowledge the value of discussions with Professor 
R. E. PEERLS. I am also indebted to Mrs. I. EGoRrOv for help with numerical 
calculations and to members of the staff of the Department of Mathematical 
Physics, University of Birmingham for help with the English text of this paper. 


RIASSUNTO (*) 


LS 5; . . . . 
E stato fatto l’esame delle sezioni d’urto per le reazioni 
20 + nt 


Yai Ol | su 


in prossimità della soglia, Si presenta il confronto coll’esperienza dei risultati di tale 
esame e si dimostra che il rapporto negativo/positivo è costante nel campo d’energia 
fotonica (153-193) MeV ed uguale a 1.4 (con l’approssimazione del 10%), d’accordo 
coi risultati della teoria dei campi. Il disaccordo dei precedenti risultati col rapporto 
di Panofsky e gli spostamenti di fase dello scattering delle onde S si discute dal punto 


di vista dell’ipotesi del campo dei mesoni 7°. Si propone una prova sperimentale diretta 
per la verifica di questa ipotesi. 


(*) Traduzione a cura della Redazione. 
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The Interaction between Equilibrium Defects 
in the Alkali Halides: 
the « First Excited State» Binding Energies 
of the Impurity Complex and of the Vacancy Pair. 


VE ae OSL (2) 
Istituto di Fisica dell Universita - Palermo 
and 
G. ATROLDI 
Laboratori CISH - Milano 


(ricevuto il 1° Marzo 1958) 


Summary. — The experimentally observed presence in NaCl and KCl 
crystals of comparable concentrations of complexes between divalent 
impurity cations and positive-ion vacancies in the ground and first excited 
states, is Shown to be in agreement with the results yielded by the Born- 
Mayer theory of ionic crystals. Theoretical estimates for the binding 
energy of the first excited state of the Sr?+ impurity complex and of the 
vacancy pair in NaCl and KCl erystals are presented and discussed. 


1. — Introduction. 


Divalent impurity cations and cation vacancies, as well as anion and cation 
vacancies, attract each other in an alkali halide lattice as they bear opposite 
effective charges, and one may expect a priori that pairs of these oppositely 
charged point imperfections will be present in the crystal lattice in various 
states of association, corresponding to various distances of the two interacting 
defects. While the ground state of these pairs of point defects in the alkali 
halide crystals has been subjected in recent years to several theoretical in- 


(*) Present address: H. H. Wills Physics Laboratory, University of Bristol, England. 
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vestigations (**), only preliminary results have yet been reported for the ex- 
cited states (4). The point is that experimental data have been available for 
some time on the ground state of these associations, but no experimental 
method had been found till recently for a comparative study of the various 
states of association between point defects. Only recently the use of electron 
spin resonance techniques has allowed the determination in NaCl and KCl 
crystals of the relative concentrations of the lower association states of com- 
plexes between aliovalent paramagnetic impurity cations and cation va- 
cancies (*). The principle of the method ($) is that a vacancy present in the 
neighbourhood of a paramagnetic ion reduces the symmetry of the crystalline 
electric field acting on the ion in a way which depends on the relative position 
of the two defects in the crystal lattice, and the dependence of the absorption 
spectrum on the orientation of the crystal relative to the static magnetic field 
is strictly tied to the symmetry of the crystalline electric field. 

In this paper we estimate theoretically the binding energy of the first ex- 
cited state of the impurity complex between Sr?* ions and positive-ion va- 
cancies in NaCl and KCl crystals: this binding energy turns out to be com- 
parable with the binding energy of the ground state, suggesting the model 
of a spread out impurity complex. We evaluate also the binding energy of 
the first excited state of the vacancy pair in the same crystals, which is found 
to be smaller than the binding energy of the ground state by a few tenths 
of an eV. 


2. — The first excited state of the impurity complex. 4] 


While the interaction energy of two charged point defects in the alkali 
halide lattice, in the limit of large relative distance, is equal to the Coulombic 
interaction energy of two point charges in a continuous medium characterized 
by the static dielectric constant of the crystal, the Born-Mayer model of ionic 
solids yields a value rather less than Coulombic for the binding energy of the 
ground state of the impurity complex (?) and of the vacancy pair (*) in NaCl 


R. Reitz and J. L. GAMMEL: Journ. Chem. Phys., 19, 894 (1951). 
. BAssaNI and F. G. Fumi: Nuovo Cimento, 11, 274 (1954). 
M. P. Tost and F. G. Fumi: Nuovo Cimento, 7, 95 (1958). 
. BassaNI and F. G. Fumi: Suppl. Nuovo Cimento, 1, 114 (1955). 
G. D. Watkins and R. M. WaLKkER: Bull. Amer. Phys. Soc., Ser. II, 1, 324 
ce and private communications. The same method has been applied to NaF crystals 
by W. Hayes (unpublished), and to MgO crystals by J. E. WeRTZ and P. Avzins: 
Phys. Rev., 106, 484 (1957). 

(6) E. E. ScHNEIDER and J. E. CAFFYN: Report of the Conference on Defects in Crys- 
talline Solids held at Bristol in July 1954 (London, 1955), p. 74. 


33 - Il Nuovo Cimento. 
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and KCl crystals. In fact, the polarization. of the material between the two 
defects increases their interaction energy, while the polarization of the material 
outside the defects decreases it (7). On these qualitative grounds one may 
expect that the binding energy of the first excited state of the impurity complex 
in these crystals will be nearly equal to the Coulombic value, i.e. comparable 

with the binding energy of the ground state 


re as calculated by BASSANI and FUMI (?). Indeed, 

7 4 4 i è the polarization of the negative ion sitting 
: À P 5 7 between the impurity and the position where 
the cation vacancy is to be created (Fig. 1) is 

FR È very large owing to the large polarizability of 
> the ion itself and to the strong electric field due 
siberiana i to the impurity ion, and it is very effective in 
- n fi n s lowering the energy to create the vacancy 


ME nce Le 
ig. 1. — The first excited state aa NE: E ec : 
eo The binding energy of the first excited state 


cation vacancy complex in the of the impurity complex can be actually evalu- 
NaCl lattice. The approxima- ated using the Born-Mayer theory of ionic 
tion adopted for the distortion crystals, as first applied to the study of point 
of the lattice around the comp- imperfections in the alkali halides by Mort 
lex is shown. à ae à È 

and LITTLETON (5). This binding energy is 

given by the difference between the energy 

to create a positive-ion vacancy in the perfect crystal lattice and the 
energy to create it in the position of fourth neighbour of a divalent impurity 
cation. To account for the polarization of the lattice around the vacancy, 
one evaluates the energy to create a vacancy as the negative of the 
average of the potential energies in the position of the ion to be removed before 
and after it has been removed. One needs therefore the distortion of the 
lattice around the divalent impurity cation and around the impurity complex 
in the first excited state (Fig. 1). The displacement &r, of the ions which are 
nearest neighbours of the isolated impurity cation is determined by imposing 
the equilibrium condition to the electrostatic and repulsive forces acting on 
each ion (*!). The displacements €7, and & 7, of the ions which are nearest 
neighbours respectively of the impurity cation and of the cation vacancy tied 
in the complex are determined in the same way by considering only the effect 
of the nearest defect but neglecting the long range polarization since the 
complex is neutral (*!). Finally, the displacement nn of the negative ion 
sitting between the impurity and the vacancy is determined by imposing 
the equilibrium condition to the electrostatic and repulsive forces acting on 


(7) F. G. Fumi and M. P. Tost: Faraday Soc. Disc., 28, 92 (1957). 
(6) N. F. Morr and M. J. Lirrreron: Trans. Faraday Soc., 34, 485 (1938). 
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it, due to the two defects and to the polarization induced by the defects on 
the nearest ions. When the lattice distortion around the impurity and around 
the complex in the first excited state is known, the energy to create a 
cation vacancy in the position of fourth neighbour of the impurity can be 
easily evaluated following the procedure already applied in the analogous 
problem of the ground state for the impurity complex (*?). The numerical 
results obtained for Sr?* ions in NaCl and KCl crystals are collected in 
Table I. For the constants entering the calculation we have used the numerical 
values quoted in ref. (?), Table I, except that we adopt for the electronic po- 
larizabilities the values reported by TESSMAN, KAHN and SHOCKLEY (°) for 
ions in crystal lattices (1°). The repulsive form given by BORN and MAYER (11?) 
has been used. 


TABLE I. — Binding energy of the Sr?* impurity complex in the first excited state in NaCl 
and KCl erystals and related ionic displacements. 


£ | E | 5) | n Binding energy (eV) 


| NaCl 0.033 0.042 | 0.086 | 0.090 | 4.78 — 4.37 = 0.41 
Cl 0.099 0.109 0.084 | 0.164 | 4.44—3.95 = 0.49 


The approximation adopted for the distortion of the lattice around the comp- 
lex is probably the principal source of error in our calculation. However, it has 
been verified that reasonable changes of the displacements around the quoted 
values have a small effect on the binding energy: this is true also for the dis- 
placement n", which may be affected by a fairly large error owing also to the 
uncertainties of the repulsive interaction law at short distances (7). The dipole- 
dipole interaction has been found to give a substantial (negative) contribution 
to the binding energy of the first excited state of the complex, contrary to 
the case of the ground state where it can be neglected (1): this contribution 
is substantially greater in NaCl than in KCl. 

The binding energies in NaCl and KCl crystals between a Sr?* ion and a 
cation vacancy in the first excited-state (0.41 eV in NaCl; 0.49 eV in KCl) 


(9) J. R. Tessman, A. H. KAHN and W. SHOCKLEY: Phys. Rev., 92, 890 (1953). 

(10) This entails some numerical changes in the displacements (cfr. Table I and 
ref. (2), Tables III and IV) and in the energies to create cation vacancies (for 
the energy to create cation vacancies in the perfect crystal lattice, see ref. (°), 
Sect. 3). However, the values reported by Bassani and Fumi (*) for the binding energy 
of the ground state of the impurity complex remain practically unchanged (see ref. (3), 
footnote (15)). 

(11) M. Born and J. MAYER: Zeits. f. Phys., 75, 1 (1932); J. MAYER and H. HELM- 
HOLZ: Zeits. f. Phys., 75, 19 (1932). 
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are found to be comparable with the binding energies in the ground state 
(0.43 eV in NaCl; 0.38 eV in KCI) (ref. (*), footnote (16)), the ground state 
being slightly more strongly bound in NaCl and less strongly bound in KCl. 
These results are in good qualitative agreement with the experiments by WAT- 
KINS and WALKER (5), who observe by paramagnetic resonance techniques 
an inversion in the relative binding energy of the ground and first excited 
state for the Mn2+ impurity complex in going from NaCl to KCl crystals: 
the ground state is more strongly bound in NaCl, and the difference in binding 
energy between the two states is quite small in both salts, of the order of a 
few hundredths of an eV (2). Moreover, our model of a spread out impurity 
complex can explain qualitatively ASCHNER’s observations which indicate that 
the associated cation vacancies are as effective as the free vacancies in aiding 
the self-diffusion of Na+ ions in Cd?+-doped NaCl (13:14). 

It is also worth noting that the value yielded by our lattice calculation for 
the binding energy of the first excited state of the Sr?* impurity complex is 
comparable with the Coulombic value (0.46 eV in NaCl; 0.49 eV in KCl), so 
that the approximation adopted by LIDIARD (!) in the analysis of the ionic 
conductivity of the alkali halides doped with divalent impurity cations ap- 
pears justified. 


3. — The first excited state of the vacancy pair. 


To calculate the binding energy of the first excited state of the vacancy 
pair, one has to evaluate the energy to extract a positive ion from the position 
of third neighbour of a negative-ion vacancy, and to subtract it from the 
energy to extract a positive-ion from the perfect crystal lattice. The 
determination of the distortion of the lattice around the vacaney pair 
in its first excited state is straightforward: the outward displacements &r, 
and &,r, of the nearest neighbours of the positive- and negative-ion vacancy 
are evaluated as the displacements of the nearest neighbours of the isolated 
vacancies, neglecting the long range polarization. The numerical results 
are collected in Table IT. It can be noted that the Born-Mayer model yields 
for the binding energy of the first excited state of the vacancy pair a value 


OSES TE CAAT ANNAN NAAN NTT SANS Sert 
è Wa | 


(2) It may also be expected that in NaF, where the polarizability of the anion 
is considerably smaller than in NaCl, the ground state of association should be more 
strongly bound, as observed by HAYEs (private communication). 

(18) J. F. ASCHNER: Self-Diffusion in Sodium and Potassium Ohloride, Ph. D. Thesis 
(University of Illinois, 1954). 

(4) A. B. LiprarD: Handb. d. Phys., 20, 246 (1957), and private communication. 

(5) A. B. Liprarp: Phys. Rev., 94, 29 (1954). 
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TABLE II. — Binding energy of the first excited state of the vacancy pair in NaCl and 
KCl erystals and related ionic displacements. 


| 


Ja) | S | Sa | Binding energy (eV) | 
| | | 

NaCl 0.104 0.086 0.129 | 4.78 — 4.50 = 0.28 
KCl 0.076 0.084 | 0.095 4.44 — 4.06 = 0.38 | 


(*) The values of the quantity ESA which is the displacement of the nearest neighbours 
of an isolated anion vacaney, are given by Tosi and Fumi (ref. 3, Table I). 


substantially smaller than for the binding energy of the ground state (0.60 eV 
in NaCl; 0.72 eV in KCI) (3). Moreover, the value is still definitely smaller than 
the Coulombic value (0.53 eV in NaCl; 0.57 eV in KCl), as can be easily under- 
stood on the basis of the qualitative argument given at the beginning of Sect. 2. 


The authors are indebted to Professor F. G. FUMI, who suggested the problem, 
for helpful discussions and continued encouragement. Very informative cor- 
respondence with Dr. G. D. WATKINS is acknowledged. 


RIASSUNTO 


Si mostra che la presenza in cristalli di NaCl e KC) di concentrazioni confrontabili 
di complessi fra cationi d’impurità bivalenti e posti vacanti da ione positivo negli stati 
fondamentale e primo eccitato, osservata sperimentalmente, è in accordo con i risultati 
forniti dalla teoria di Born-Mayer dei cristalli ionici. Sono presentati e discussi risultati 
teorici per l’energia di legame del primo stato eccitato del complesso di impurità per 
ioni Sr?+ e della coppia di posti vacanti in cristalli di NaCl e KCl. 
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Charge Conservation in the Bethe-Salpeter Equation. 


G. R. ALLCOCK 


Department of Theoretical Physics - University of Liverpool 


D. J. HooTON 


Department of Applied Mathematics - University of Liverpool 


(ricevuto il 5 Marzo 1958) 


Summary. — The scalar product between two bound states is expressed 
in terms of a conserved charge, and reduced to a standard covariant form. 


1. — Introduction. 


In a previous paper (!), hereafter referred to as I, one of us has derived 
a general formula expressing the scalar product between two bound states 
in terms of their Bethe-Salpeter wave functions. MANDELSTAM (?) has derived 
an alternative form for the scalar product, applicable to states where there is 
a non-zero conserved charge (e.g. electric charge, nucleon charge, ete.), and 
this form has been used by various authors (3). 

The two forms for the scalar product are quite different in appearance, 
though of course formal considerations indicate that they should be equivalent. 
In the present paper we treat a system in which the conserved charge is carried 
by one or more particles, and trace out explicitly the equivalence between 
the scalar product derived from charge quantization and that given in I. 
We also indicate that when the Bethe-Salpeter kernels and one-particle pro- 
pagators are approximated by dropping out various graphs the equivalence 


(4) G. R. ALLcocK: Phys. Rev., 104, 1799 (1956). 

(?) S. MANDELSTAM: Proc. Roy. Soc. (London), A 233, 248 (1955). 

(?) K. NISHIJIMA: Progr. Theor. Phys., 13, 305 (1955); A. KLEIN and C. ZEMACH: 
Phys. Rev., 108, 126 (1957); F. L. Scarr and H. UMEZAWA: Phys. Rev., 109, 1848 (1958). 
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persists, provided a proper correspondence between graphs is maintained 
between the two methods. 


2. — Matrix elements of the charge. 


We consider two bound states <b| and |a), described by Bethe-Salpeter 
wave functions y, and y,, as in I. We suppose, for the sake of the argument, 
that the conserved current j” is carried by a Dirac field y, so that it takes 
the form py“y. The current matrix element is then 


(1) <b |9"(x)|a> = <b|y(x)y“y(x) |a>d . 
MANDELSTAM (?) has shown that it can be written covariantly in the form 
(2) COG) == tee eyes 


Here we imply integrations over the space-time arguments of the Bethe- 
Salpeter wave functions % and y,. These space-time arguments appear also 
in the graphic expression R(x), which arises as follows. 

To express the current matrix element (1) in the form (2) one augments 
the operator arguments of the propagator G descriptive of the bound system, 
by putting in the current operator %(x)y"y(x). Those graphs belonging to 
this new propagator and contributing to (1) may be enumerated one-to-one 
by inserting a vertex y" in turn into each w-line of the graphs for G. Since G 
has the structure expressed by the equation 


(3) G= G+ GAG, 

the propagator with current has the structure 

(4) GR*(x)G, 

where R(x) in turn can be written as a sum of two terms, 

(5) (ERE) Gr) KE 
Here K"(x) and G*(x) are the expressions arising from vertex insertions in 
the y lines of K and G, respectively. 

Inserting (5) into (2) and integrating over all space we obtain the matrix 


element of the total charge, 


(6) iQ |a) = deg (Ko(w) + KGSw)K) xa - 
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Mandelstam’s (2) scalar product is obtained by exploiting the usual rela- 
tionship between the charge operator Q(t) and its integer eigenvalue n. Thus 
for the scalar product one should have 


‘I 
(7) b|a) == <b|Q(t)|a>, 


n 


where <b|Q(t)|a) is given by (6). 


3. — Simplification of the charge matrix elements. 


If we analyse % and y, into momentum eigenstates, the charge matrix 
element (6) appears as a sum of contributions, one from each momentum 
eigenvalue. For the space integration in (6) gives momentum conservation. 
It follows that when the states |a) and |b> have the same mass eigenvalue M 
the time dependence from y, in (6) is completely compensated by that in %, 
so that (6) becomes explicitly independent of the time t, irrespective of the 
formal theorem of charge conservation. Naturally our interest in the present 
work is confined to this case, since with unequal masses <b|a@> vanishes by 
space-time translational invariance. 

In view of this explicit time independence it is legitimate to replace (6) 
by its average over a long time interval (7, — 7); thus 

7 
(8) IQ |a> = sp farm) + KOM} 70. 
n° 
The above expression corresponds in some sense to a charge measurement in 
which the energy of the system is completely undisturbed, and leads to con- 
siderable simplifications. For the four-dimensional integrations on æ in (8) 
allow us to apply Ward’s equation 


(9) SP) Sp) = — OS,(p)/Op, 5 
or 
(10) [Sete — £)d4ey" SG — y) = = (af — VINS (T1 — vi). 


In the matrix element (8) we do not actually encounter the expression (10), 
but rather the expression 


r 
ia 2 
(11) oT iS,(@, — ®) dtay" iS p(a —Y,). 


aft 
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Here #, and y, are the end points of a section of a y-line into which a vertex y" 
has been inserted. To see how to deal with (11), we must consider the in- 
tegrations over the space-time co-ordinates other than æ in (8). These integ- 
rations extend over the whole of space-time, but for any pair of states |b) 
and |a) we will be able to find some large but finite time interval 27 surround- 
ing #°, outside which the integrand is effectively zero. This being so, the 
points x, and y, with which we have to deal in (11) will be separated in time 
by an amount of order 27, or less. If then we take 7 > 7, the time integ- 
ration in (11) usually extends far to each side of the interval (x°, y?), and the 
contribution of the pair of points v,, y, can then be calculated by using (10). 
Therefore we may with sufficient accuracy replace the expression (11) by 


ul 
(12) 27 
0, x, and y, not both within (7, — 7). 


(ay — ya — Vv). x, and y, both within (7, — T), 


The substitutions (12) are correct apart from edge effects arising when x, or y, 
lie within a time range 7 of the limits --7. The edge contributions appear 
multiplied by a factor 1/27 from (11), and therefore vanish as t/T as T — co. 
Hence, if 7/7 is sufficiently large, the integrals on x in (8) may be performed 
by replacing in turn each y-line factor iS,(x, — y;) in K or & by (12). Thus, 
neglecting edge effects, each vertex insertion brings in a factor (1/27): (a! — y#) 
if x) and y, lie within (7, —T), and a factor zero elsewhere. 

Again bearing in mind the localization within the range t, we see that to 
within edge contributions the factors x — y sum to zero for vertex insertion 
in closed loops, while for any line running right through a graph they sum 
to the difference of the end co-ordinates of the line. 

Evidently for the charge matrix element (8) we now have the expression 


ate 1 a 
(LI) cu il tr fat Got. More == URW lees A aaa cca BP i 
SO) (UE Une) 


Here we understand that X and G, are decomposed graphically so that the 
charge-carrying lines can be seen. The co-ordinates x, and y, denote the end 
points of a typical charge-carrying line, and > denotes a summation over all 
such lines. Integrations over the whole of space-time are understood for the 
co-ordinates other than «x, and y,. 

At this stage it is useful to introduce centre co-ordinates X#, Y" and re- 
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lative co-ordinates rf, sf thus: 


v 


| Dé = > CAE « yey = > CAE > a; = il , 
(14) | L L i 


o — al qu Ae Re 7 
ri Xe, CAI — a 


Here the «; are arbitrary real numbers introduced only as an aid in the cal- 
culations. We now substitute 


(15) at — y) = Xo+ ri — (¥°+ st), 
and apply the Bethe-Salpeter equations 7 = GKy, = y~KG, to the second 
term on the R.H.S. of (13). Thus we have, for example, 
a 

(16) rt fan Gu. COSTOSO cent «WES Patios I) 

7 = 7% {al % ...) + edge terms}. 
The edge terms above arise only when «} lies within a range t of +7, and 
they are multiplied by the factor r°, which is confined to within the range 
(t, — t) by the wave function y,. Evidently these edge terms can be neglected, 
because of the factor 1/27 in (13). A similar application of the Bethe-Salpeter 
equation to the term X° | G,K 4% is useless, giving non-vanishing edge errors 
because X° is of order +7 in the edge regions. Substituting (15) in the 
first term on the R.H.S. of (13), and (16) in the second term, the contributions 


of the relative co-ordinates r and s cancel between the first and second terms 
in the limit 7 — co, and we are left with 


di; TI, 
| ue: 
(LTD Q ay = > li sp fate fanta. IRR VO) 
—T —T 


Et es Ve OOO 0 TUE e 
RE NOR OA oo 


The contributions to (17) come from the region — 7 < X°— Yon 7, This 
being so, we may with negligible edge error replace 


Ti (co) (co) 
fas | dif dr, 350 
7 — oo) 


(col 


148 
fata ay, 
FA) — 0 
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mit 


where 


(18) ddr) dir, dir, dre 


i 


The y integrations can be treated similarly. When this is done the summation 
over the charge carrying lines becomes trivial, since the factors X° — Y° are 
the same for all lines. If there are n+m forward directed charge-carrying 
lines and m backward directed charge-carrying lines the summation brings 
simply a factor n. Since the formal theorem of charge quantization ascribes 
to <b|Q|a> the value n<b|a>, we obtain now for the scalar product (n # 0), 


ni vi 
ANA = 
(LO) echi aa Tuna sp [ex faro. RO AO 
=P =f 


Eee) 
Gn ee E AA 


Integrations over the elements dV (eq. (18)) accompany the X and Y in- 
tegrations. 

In a theory in which all propagators and wave functions are defined with 
respect to renormalized Heisenberg operators Lo, eq. (19) still holds, for 
<b|Q|a> then should have the value Z,'n<b|a), while vertex insertion brings 
a factor (1/27)Z,, "(x — yf) because of the renormalization vertices — 1(1 — Z,)- 
«(0 — m). 

The integrations on the centre co-ordinates in the scalar product (19) must 
now be performed. Since X° — Y° is confined to a range 27 < 27 by the 

Chi a di (co) 

wave functions y, and %, we may replace I axe faye by fax fa — Yo) 
with vanishing edge error. If we now perform all the integrations other than 
that on X°, the translational invariance ensures momentum conservation, so 
that the frequencies in y, annul those in y,, giving an integrand which is 
independent of the integration variable X° except in the edge region. The 
integration on X° gives therefore merely a factor 27 to cancel the factor 1/27 
ina (10) 

To express this feature explicitly we introduce Fourier transforms for G, 
K, y and x, as in I, thus 


D 
2 
8 
SS 
I 
© 
Ss 
ay 
jon 
i 
es} 
R 
S 
S 
% 
© 
bd 
iS 
| 
d. 
cS 
[= 
| 
ke 


Yale Bu.) = eur OPP) (P: Mere exp (iP x. 
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Using the above Fourier transforms, and suppressing the relative co-ordinates, 
the scalar product (19) now becomes 


0G(P) 
(21) (bla Ca) E 


sante, t K(P 
apo oP.) GP, Le 


es > 1 fèK() KO) wr). 


The factor 0K(P)/0P, arises because we deal with the Fourier transform of 
ON AE Ca EE NT IS «== 001 (Ds cca eee ese) CD 
The factor 0G,(P)/0P, arises in the same way. 

An expression very similar to (21) was obtained by a different method in I. 
In the following section we shall show that the two expressions are, in fact, 
explicitly equivalent. We may close the present section by noting that our 
proof of the equality of the scalar product (7) and the expression (21) is valid 
even when one approximates G, and K by using subsets of the full series of 
graphic contributions. It entails only that one take just those graphs for 
the current density j“ which arise from vertex insertion in the graphs actually 
used to represent G, and Æ. For this reason computations made on the basis 
of charge quantization will necessarily agree exactly with any based on (21), 
provided the Bethe-Salpeter equations y= G,Ky are solved exactly with the 
approximated G, and K (*). 


4. — Lorentz invariance of the scalar product. 


The method used in I to determine the scalar product rested essentially 
on finding an expansion of the propagator G(p; ...r...; ...8...) in the neigh- 
bourhood of the singularity at p? = M?, M denoting the mass of the bound 
system. The value of G(p) at a point p just off the mass hyperboloid was de- 
fined by eq. (21) of that paper, which we shall quote as (1.21). The Bethe- 
Salpeter equation for G appeared in the form 


(1.22) g(p) = (p? — M?)G(p) = (p? — M?)G,(p) + GP) K(p)g(p) . 


In I the variation of the regular function g(p) across the mass hyperboloid 
was investigated by writing p“= P"(1+e), where P?= M?. We proceed here 
more generally, by considering an arbitrary displacement dp”, so that on dif- 
ferentiating (I.22) above we replace (I.25) by the four equations 


(22)  {1— G.(P)K(P)}d9(P)/dPt = 2P,G(P) + (0{ @(P) K(P)}/0P*) g(P) . 


(*) KLEIN and ZEMACH (*) have examined the form of the current density (2) in 
the lowest ladder approximation. We wish to point out that their scalar product (95), 
obtained for this case, contains one too many integrations on the centre co-ordinates. 
The dependence on the relative co-ordinates expressed by eq. (96) is correct, however. 
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NI 


As in I, we express g(P; ...r...; ...8...) as the linear combination 
(23) OP) SUD Tr.) Agg(P)Og(Pj 8.) 


where w, and v; are complete sets of linearly independent solutions of the 
Bethe-Salpeter equations 


Then, as noted in I, eqs. (22) serve to determine the matrix of coefficients 1,,. 
For the operator 1 — G,K has zero determinant, so that (22) can only be solved 
for 0g(P)/0P" if pre-multiplication of its R.H.S. by the null vectors 0,(P) K(P) 
gives zero. Hence, using (23), we obtain four sets of equations, replacing (1.27). 
They are 


(25) N gl) = 0,(P)(K(P){ 0G,(P)/AP*“\K(P) + 0K(P)/OP") u,(P) . 


Compatibility of these four sets of equations follows from the fact that 


N,xp(P) is invariably parallel to P,, as we shall presently show, 


op 


(26) N,sg(P)=P,Cag(P), 


u “xp 


so that (24) reduce to a single scalar equation 
(27) cpl) Ag (P) + 20, = 0. 


The matrix equation (27) provides a convenient specification of the coefficients 
2,, at the pole of the propagator, though it remains to be seen whether 4 
actually fulfils the requirement that the quantization metric is positive de- 
finite. 

To prove (26) we invoke the Lorentz covariance of the theory. This asserts 
that within any given frame of reference we may attribute to each infinite- 
simal Lorentz transformation of P“, P“ > P“+ SP, a corresponding trans- 
formation w,(P; seus) du;(P; ...8...). Thus, writing down the 
Bethe-Salpeter equation for u;+dv; we have, 


(28) ug + Su; = (GK + {0(G,K)/OP* dP") (ug + 3us), 
whence 
(29) (1 — GK) du; ={0(G,K)/OP*\8P"-u,. 
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Premultiplying (29) by 1K, and using the Bethe-Salpeter equation for ©, and 
noting that SP” is any infinitesimal four-vector orthogonal to P“, we at once 
obtain (26). It will be seen that this result is also independent of the sub-sets 
of graphs taken to represent G, and K. 

Returning to the normalization integral, we require the matrix 2-1, which 
from (27) is given by 


(30) A. 


We can now evaluate (30) in various ways, the equivalence of which follows 
from the identity (26). For instance, if we take the zero’th component of (26) 
we obtain 


(31) J-1 = — No/2P, 


while if we contract (26) with the four-vector P, (assuming M?= P? = 0) 
we obtain the more obviously covariant expression 


(32) 1-1 = — N*P j2M2. 


At this point we must refer the reader to Sect. 5 of I, where the scalar product 
is expressed in terms of the matrix 4-! by eq. (1:35). Substitution of (31) 
into (1.35) verifies at once our expression (21) for the scalar product, while 
substitution of (32) gives the obviously covariant expression (1.37) below: 


È St RISE ACER) Re CURE 
(Seba | Plon] op, K(P) prati )f 42) 


The scalar product (21) derived from charge quantization is obviously 
entirely equivalent to the general covariant form (1.37), by virtue of the 
identity (26). 


RIASSUNTO (*) 


Il prodotto scalare tra due stati legati è espresso in termini di una carica conser- 
vata e ridotto in forma covariante convenzionale. 


(*) Traduzione a cura della Redazione. 
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(ricevuto il 13 Marzo 1958) 


Summary. — An interpretation of the data on the capture of 2 -particles 
by nuclei of photoemulsion has been attempted on the basis of the inde- 
pendent particle model. In particular it has been found that the prob- 
ability of the A-particle produced in the capture to remain trapped in 
the nucleus, and hence the average size of the Y capture star, strongly 
depends on the mean free path of A-particles in nuclear matter. In order 
to justify the experimental data the A-N elastic cross-section at energies 
around 50 MeV should be about 15 mb. 


1. — Introduction. 


Various calculations based on the independent particle model and with 
the Goldberger method, performed either analytically or with the help of the 
Monte-Carlo method, have successfully accounted for the mechanism of nuclear 
disintegrations induced by fast protons (1) and by the capture of m- and K- 
mesons (2%). Similar calculations have also been found useful for deriving 
information on the interaction properties of «strange particles » with nu- 
cleons (*), from the analysis of the more easily available experimental data 
on complex nuclei. 


(1) G. BERNARDINI, E. T. BoortH and 8. J. LINDENBAUM: Phys. Rev., 88, 669 (1951). 

(2) V. DE SABBATA, E. MANARESI and G. Puppi: Nuovo Cimento, 10, 1704 (1953). 

(3) F. C. GILBERT, C. E. VioLet and R. S. Wire: Phys. Rev., 107, 228 (1957). 
R. H. Gapps: Phys. Rev., 107, 239 (1957). 

(4) See for instance the papers of the Gottingen group on K* and K° nuclear 
interactions as well as of the Bologna, Bristol-Dublin and Padua groups on K* inter- 
action: Nuovo Cimento and Zeits. f. Naturfor. (1956-1957). 
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This encouraged us to apply the same method for the nuclear capture of 
=> hyperons. The purpose was to give an account of the experimental data 
on this type of events and, more specifically, to derive some information 
about the nuclear interaction of A particles of which the 27 capture process 
constitutes a convenient source. 


2. — Disintegration induced by the &% capture. 


The capture of a 27 particle by a nucleus will occur from a state of orbital 
angular momentum which depends on the cross-section of the X-nucleon inter- 
action. For the present calculation we have assumed that the cross-section 
for the process YX+N +—A+#N is rather small (~ } geometric, as indicated 
by the work of FERRARI and TONIN (5) and of ALLES (*) so that the X is able 
to penetrate inside the nucleus before being actually absorbed. Consequently we 
have taken the points of capture to be uniformly distributed inside the nucleus. 
It can be seen however that a change of this distribution, for instance cor- 
responding to a Coulomb potential oscillator or to capture occuring only at 
the periphery plays little role in the final results of the calculations. 

The X-particle may be captured following the two reactions: 


a) Xx +p->A+n4+81 MeV 


b) x +p=>X°+n+ 8 MeV 


plus, eventually, reactions involving two nucleons. Since in the latter case 
both the nucleons have to be emitted with energies higher than the maximum- 
Fermi energy, these reactions are likely to be strongly depressed by the Pauli 
exclusion principle. For this reason, and also in analogy with K-meson 
capture (7), two nucleon reactions have been considered rare and therefore 
neglected. Of the two reactions involving a single nucleon, reaction b) is found 
to be forbidden in about 90% of the cases by the Pauli principle because of 
the low Q-value. Further, out of the two subsequent reactions: 


€) ZÒA+y 
and 
d) > NS Ae 


(5) F. FERRARI and M. Tonin: Proc. of Padua-Venice Conference (1957). 

(6) W. ALLES: private communication. 

(7) According to different groups the frequency of two-nucleon K~ captures lies 
between 3 and 10 percent. (See Padua-Venice Conference report). 
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of the Y° created inside the nucleus, the second one is likely to be dominant. 
Since the energy spectrum of the A-particles emitted from the reaction d) 
will be practically the same as in the reaction a), we may consider only the 
latter one as responsible for the development of the A cascade inside the 
nucleus. 

The spectrum of A-particles emitted in reaction a) has been computed 
with essentially the same method used by Capps (5) for obtaining the spectrum 
of X-particles in the reaction 


Re NERE 


The depth of potential well, V, of the A-particle has been taken to be 17 MeV, 
as resulting from the formula 


V= 0(3U, + U,)/4, 


where o is the nuclear density and U,, are the volume integrals of the A-N 
interaction potential for parallel and anti-parallel spins. For these we have 
taken the values 84-10-%° and 403-10-3° MeV 
em’ calculated by DALITZ (*) from the binding 
energies of light hyperfragments. The energy 
spectrum of the A-particles inside the nucleus 
is given by the curve of Fig. 1. 

From the point of capture we have then to 
trace both the nucleon and the A-particle up to Se a tO 
the nuclear surface, where they will either escape A (in) Fe 
or be trapped if the kinetic energy is lower than pis. 1.- Energy spectrum at 
the corresponding potential step; this has been production of A-particles inside 
taken as energy independent. the nucleus. 

Due to its great complexity this part of the 
calculation has been done with the help of the Monte-Carlo method. For the 
cascade induced by the neutron, the values of the various parameters (o,,,, 
Oi ete.) have been taken to be the same as in reference (?); for the A cascade 
the A-N scattering in the center of mass system has been assumed to be iso- 
tropic and energy independent, the total cross-section on free nucleons, 
O,x, being left as a free parameter. 

Four sets of calculation have been performed with different values of o, x, 
namely 0, 15, 30 and 60 mb. The cross-sections on protons and neutrons have 
been taken to be equal as expected from the fact that the isotopic spin of 


N 


a 


m 


d 


(3) BR. H. Capps: see ref. (3). 
(*) R. H. Darirz: Rep. on Progr. in Phys., 20, 252 (1957). 


39 - Il Nuovo Cimento. 
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the A is zero. The calculation has been done for a nucleus with A = 100. 

After the nuclear cascade, the nucleus will be generally left in an excited 
state and nuclear particles will evaporate. In addition also the A-particle 
previously trapped may evaporate; we have however reasons to believe that. 
this will occur rather infrequently. Indeed, if the A-particle has a non-vanishing 
cross-section with nucleons, for instance around 10 mb, after being trapped 
it will collide with a number of nucleons in a time of the order of 1072! s. The 
effect of each collision will be, as a consequence of the Pauli exclusion prin- 
ciple acting only on the nucleon, to reduce the A energy; therefore in a very 
short time, the A-particle will almost reach its ground state. Its subsequent 
re-acceleration and evaporation should then be even less probable than the 
evaporation of a given neutron of the nucleus: a process which at the excitation 
energies here concerned seems to take place with a probability of only about 5%. 
This reasoning has led us to neglect the evaporation of the A-particle and 
only to treat that of nucleons; this has been done with the same method as 
Inarererence:(?). 

After a time much longer than the evaporation process the trapped A-part- 
icle will undergo decay. In a heavy nucleus mesonic decay should be practi- 
cally absent and the dominant reactions will be 


e) AX nen 
and 
f) ARDEA Die 


Experimental data on the relative frequencies of reactions e) and f) are still 
searce for light nuclei (7°) and totally absent for the heavy ones. In our cal- 
culations we have assumed the two reactions to occur with equal probability; 
for o,.=15 mb additional computations have been made for the cases of 
decay « stimulated » exclusively by protons, or by neutrons. 

The two nucleons emitted from the reactions e) and f) have been followed 
through the nucleus and the emission of particles due to this second cascade 
and in the consequent evaporation has been computed. 

The calculations have been done mostly graphically with the help of a 
device which permits the construction of the model of a three dimensional 
nuclear cascade. For the cascades initiated by a nucleon as well as for com- 
puting the emission of evaporated particles, we have had the opportunity of 
utilizing a number of partial data of the Monte-Carlo calculation performed by 


(1) M. BaLDO-CEOLIN, C. DizwortTx, W. F. Fry, W. D. B. GREENING, H. Hoziva 
S. LIMENTANI and A. E. SICHIROLLO: preprint and also reported at the Padua-Venice 
Conference (1957). 
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DE SABBAT A, MANARESI and PupPI for the 7--meson capture. The number 
of «events » on which our calculation is based is as follows: 435 for o,, = 0, 
217 for o\y= 15 mb and 109 for each of the cases of o,y= 30 mb and 60 mb. 

In Fig. 2 the points linked by 
LT lines give the calculated fraction 
i pat of capture stars with a given 


Fig. 2.— Prong distribution of =~ Fig. 3. — Visible energy release of 
capture stars for different values x capture stars for different values 
of the A-N cross-section. OLO 


number of prongs. In Fig. 3 is shown the distribution of the total visible 
energy plus the binding energy of the capture stars with at least one prong. 

The fraction of A-particles trapped in the nucleus for the four cases of 
Fig. 2 and 3 are found to be 0.0, 0.20, 0.42 and 0.58. 


3. — Comparison with experiment. 


The results of this calculation have been compared with the data of these 
four nuclear emulsion groups: 
University of California (Livermore group) (11); 
Max-Planck Institut für Physik, Gottingen (2); 
University of Wisconsin (1); 
University of California, Berkeley (Barkas group) (1). 
The first two groups have attempted a direct determination of the per- 
centage of © -capture stars in which no ionizing prongs are emitted (1). Both 


(11) See ref. (3). 

(12) W. ALLes, N. N. Biswas, M. CECCARELLI and J. CRUSSARD: Proc. Padua- 
Venice Conference (1957). 

(13) W. F. Fry, J. ScanEPs, G. A. Snow, M. S. Swami and D. C. Wozp: Phys. 
Rev., 107, 257 (1957). 

(14) W. H. BARKAS: VII Rochester Conference Report (1957). 

(5) See ref. (11) and (12) and also WHITE's report at the VII Rochester Conf. (1957). 
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results, although obtained with different methods agree well estimating that 
zero prong © capture stars should be at least half of the total. The com- 
pilation of these two statistics gives the corresponding point (extreme left 
point) of Fig. 2. For the other points of this figure we have included the data 
of the Wisconsin and Berkeley groups as well as some more data of the Got- 
tingen group, of good statistical weight but reliable only for captures with at 
least one prong. The data correspond to a mixture of captures by heavy and 
light nuclei; the implications of this will be discussed in Sect. 4. For the data 
of Fig. 3 we have instead attempted to reduce the contribution of captures 
by light nuclei by eliminating as usual those events showing a proton ot less 
than 4 MeV or an «-particle of less than 9 MeV. 

Fig. 2 and 3 indicate that a complete transparency of the nucleus to the 
A-particle can be excluded because then the occurrence of a number of stars 
with several ionizing prongs or with rather high total visible energy, which 

needs to be interpretated as due to the 

AN trapping of the A, could not be explain- 
| saa teal stimulated by oi ed. On the other hand, a cross-section 
| a (dy 218 mb) | around «geometric value» (a(h/m_c)? = 
= 63 mb) or higher would not lead to 

an agreement with the observed high 

ae percentage of small stars. The curve cor- 
e ase % responding to o,, =15mb seems to fit 
Evis *E bing MeV better the experimental points. For this 

Fig. 4. — Visible energy release of X° Case, as mentioned previously, we have 
capture stars for A-decay stimulated also computed the visible energy distri- 
exclusively by protons or by neutrons. butions, shown in Fig. 4, for the decay 
of the A stimulated exclusively by a pro- 

ton or by a neutron. Since these curves differ little, the uncertainty on the 
ratio between the two reactions is not likely to be critical for our calculation. 

However, the relative frequency of the two reactions e) and f) should be 
substantially reflected in the energy spectrum of individual prongs of energy 
of more than about 40 MeV. Present statistics in this energy range is not 
sufficient for any conclusion. 


4. — Effect of captures by light nuclei. 


In the previous section the results of a calculation based on a model to be 
applicable to heavy nuclei have been compared with experimental data of 
nuclear emulsion which include also events in light nuclei. The elimination 
of these does not seem at present feasible for all classes of events; for instance 
for the zero prong stars. 
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An «-particles model is probably a better picture for a light nucleus; it 
has been found for instance adequate to account for 7--captures in C, N 
and O (!). A calculation with this model (1) is however very involved and, 
as in case of x captures, is not likely to account fully for the experimental 
facts (18); however some qualitative features of the differences between X°- 
captures by light and heavy nuclei may be discussed. 

As compared to the independent particle model in which the available 
energy is shared between two baryons, in an «-particle model this will be shared 
among five baryons so that the A will be emitted with a lower average energy (*). 
However due to the shorter nuclear path available, the probability of a A to 
loose energy by scattering will be less. These two effects will therefore play 
oppositely as far as the trapping probability of the A is concerned. Further, 
the cascading of particles in light nuclei will be less frequent so that the average 
excitation energy will be smaller. The occurence of charged prongs may how- 
ever be not depressed since the Coulomb barrier is smaller. 

These considerations indicate that the behaviour of captures by heavy 
nuclei might not be drastically different from those by light ones. For the latter 
case we may expect a lower percentage of zero prong capture stars and a 
corresponding increase of stars with low visible energy, in analogy with the 
case of 7--captures in which are involved energies of the same order as here. 


5. — Concluding remarks. 


In a recent paper HORNBOSTEL and ZORN (1°) have also obtained from the 
study of the trapping probability of A created inside the nucleus after the 
capture of a K~-meson an indication that the A-N scattering cross-section is 
significantly smaller than the geometric value. Their investigation deals with 
a rather complicated chain of reactions and with high excitation energy so 
that their results seem to be less quantitative than ours. 

Assuming that at energies here considered the scattering is essentially 
isotropic, the A-N elastic cross-section may be derived from the values of the 


(1) M. G. K. Menon, M. MurRHEAD and O. RocHar: Phil. Mag., 41, 583 (1950); 
P. Amuirasu and L. M. LEDERMAN: Nuovo Cimento, 4, 282 (1956). 

(1?) A. C. CLARK and S. N. RUDDLESDEN: Proc. Phys. Soc., A 64, 1064 (1951). 

(18) M. DEMEUR, A. HuLeux and G. VANDERHAEGHE: Nuovo Cimento, 4, 509 (1956). 

(*) Compare, for instance, the nucleon spectrum for x--capture by heavy nuclei 
of Fig. I of ref. (2), having a maximum around 45 MeV (outside the nucleus), with 
that of captures by «-particles given in Figs. 2 and 3 of ref. (7), with maximum 
around 20 MeV. 

(19) J. HornBosteL and G. T. ZORN: preprint. 
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volume integrals of the interaction potential, with the formula: 


This results from the Born approximation for a spherically symmetrie square 
potential, but is also more generally valid (2°). Using the values of U, and U, 
given in Sect. 2, o,, turns out to be 14.6 mb, in agreement with the value of 
about 15 mb, previously obtained. 


* OK OK 


We are most grateful to W. ALLES, who has worked with us at an early 
stage of this research for his contribution in discussions and suggestions. 

We would further like to thank very deeply Prof. G. Puppr and Dr. A. 
STANGHELLINI, as well as Drs. J. CRUSSARD, V. DESABBATA, K. GOTTSTEIN, 
E. MANARESI, G. QUARENI and P. WALOSCHEK for valuable advices and cri- 
ticisms. One of us (N.N.B.) is grateful to the Istituto di Fisica di Bologna 
for hospitality. 


(2) H. BetHE and R. Wirson: Phys. Rev., 83, 690 (1951). 


RIASSUNTO 


Sulla base del modello a particelle indipendenti abbiamo cercato di interpretare 
i dati sperimentali sulla cattura di particelle = dai nuclei di emulsione fotografica. 
È risultato tra l’altro che la probabilità che le particelle A prodotte nella cattura 
rimangano nel nucleo, e quindi la dimensione media delle stelle di cattura dei X-, 
dipende fortemente dal cammino libero medio delle particelle A in materia nucleare. 
Allo scopo di render conto dei dati sperimentali è opportuno assumere per la sezione 
di urto elastica A-nucleone, ad energie intorno ai 50 MeV, un valore di circa 15 mb. 
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IL NUOVO CIMENTO Vor. VIII, N. 4 16 Maggio 1958 


On the Connection of Spin with Statistics. 


N. BURGOYNE 


Institute for Theoretical Physics - Copenhagen 


(ricevuto il 20 Marzo 1958) 


Summary. — The relation between spin and statistics is proved using 
only the most general assumptions. 


The purpose of this note is to give a proof of Pauli’s well known theorem 
on the connection of spin with statistics. The argument uses recently de- 
veloped techniques of quantum field theory (*?) and establishes the theorem 
in great generality. In particular, no assumptions are made about the form 
of the field equations or interactions (*). We show that if a relativistically 
invariant field theory has the properties: 


a) No negative energy states; 
b) The metric in Hilbert space is positive definite; 


e) Distinct fields either commute or anticommute for space like se- 
parations, 


then no field can have the « wrong» connection of spin with statistics. 
Consider first a field ®, transforming under some finite dimensional irre- 
ducible representation A + S(A) of the homogeneous Lorentz group (without 


(1) A. S. WIGHTMAN: Phys. Rev., 101, 860 (1956). 
(2) D. Hatt and A. S. WiGHTMAN: Dan. Vid. Selsk, 31, no. 5 (1957). 

(3) After the present work was complete, I received a preprint of a paper on the 
same subject by G. Lijpers and B. Zumino which deals with the case of spin zero 
and spin one half fields. For the case of neutral fields their result coincides with mine, 
but for charged fields they introduce a superfluos hypothesis of « gauge invariance ». 
The reader is referred to their paper for a description of previous papers on the same 
subject. 
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inversions). Then 


U(0, A) D(})U(0, A) = D(f,), 


where 


= [Ataf Aa). fia) = f(A-20)S,,(A-) , 


and the f belong to an invariant space of test functions. fa, A} U(a, A) 
is a unitary representation of the inhomogeneous Lorentz group. U operates 
on a Hilbert space of state vectors, containing a unique vacuum 9,. We 
define D*(f)—[®(f)[*, where f is the complex conjugate of f. The metric 
is ot = ao — x" 

The twofold vacuum expectation values are 


1 (4) = (2, D (a) OF (y) Qo) ; G at) = (2, D (x) B,(y) Qo) , E=xL—Y, 
F and G can be extended to functions of a complex 4-vector 2=:&— in, 
analytic when 2? varies in the complex plane cut along the non-negative real 


axis (4). The region ¢? < 0 is within the cut plane, and in this region 


CE) 


the upper sign for integral spins (S single valued), the lower for half-integral 
spins (see Appendix for proof). 
Using these results the « wrong » commutation relations lead to 


PF (8) == G, ,(&) == 0 È? < 0 
and by analyticity 


F,+@, = 0 identically . 


Now consider the state vectors D(f)Q, and O*(f)Q,, where f(x) = f(— x). 
We compute, 


ID) + | O* AQP? = [aa [Gale — y) + Fale— y)f (y) = 0. 


(4) The argument given in (1) is for scalar fields but is easily extended to the case 
of arbitrary fields. See, for example, A. $. WiGHTMAN hectographed notes on lectures 
at Institute Henry Poincaré (1957). The same remarks apply to the result that a theory 
is uniquely determined by its vacuum expectation values. 
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Therefore D(f)Q, = O*(f)Q, = 0 for all f. Using assumption c) and the results 
of HALL and WIGHTMAN on the real analytic points we conclude that any 
vacuum expectation value containing ® vanishes identically and consequently 
that such a field is zero (24). 


The author is indebted to A. S. WIGHTMAN for his valuable comments and 
advice. He also thanks Professor NIELS BoHR for the hospitality extended 
to him at the Institute for Theoretical Physics. 


APPENDIX (5) 
G satisfies the relation, 
Suit A)S (A) (E) = Ga (AÉ) | 


by standard invariance arguments one finds that G is given uniquely as a 
sum of terms of the form 


Déc ae -Op leg nemrea G(E)., (ceo erie 


where 


Siu (A) Sa (A) Pua (ou 54 Bp) = Dale see OA as AS ; g(AË) = g(é) , 


and n is always even for S single valued, and always odd for S double valued. 
On extending G to complex values, it becomes 


Coeli, 
where g(z) = g(A,z) for 2 in the cut plane, and A, belongs to the complex 
Lorentz group (with determinant +1). In particular g(£) = g(— é) for &< 0. 


This proves the original statement. 


(5) The property proved in this appendix follows from the analysis of W. PAULI: 
Phys. Rev., 58, 716 (1940). 


RIASSUNTO (*) 


Usufruendo solo delle ipotesi più generali sì dimostra la relazione tra spin e sta- 


tistica. 


(*) Traduzione a cura della Redazione. 
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Pion Multiplicity in Antinucleon Annihilation (*). 


E. EBERLE 


Centro Siciliano di Fisica Nucleare - Catania 


(ricevuto il 26 Marzo 1958) 


Summary. — Pion multiplicity in antinucleon annihilation has been 
calculated using Dyson’s hypothesis of pion-pion resonance interaction. 
The percentage of K-mesons was also evaluated supposing that their 
volume of interaction has a radius given by the Compton wavelength 
of K-mesons. Results are in reasonable agreement with experimental 
data. 


It is well known that the application of Fermi’s statistical theory (!) to 
the problem of pion production in antinucleon annihilation gives (?) for pions 
a multiplicity = 3, 4 that is not in agreement with experimental results (8) 
n. = 4.7 40.4. 

For obtaining such a high multiplicity, with Fermi’s theory, one should 
take a radius of interaction of the order of 2 or 3 times the 7-meson Compton 
wavelength; but this hypothesis is not very reliable. 

Many efforts have been made in attempting to modify Fermi’s theory 
efficiently (*) (4); most of these were unsuccessful or not well justified. 


(*) A preliminary report on the present investigation was communicated by M. CINI 
to the International Conference of Padua-Venice (September 1957). 

(+) In a recent work of Z. KoBa and G. TAKEDA (not yet published) the authors 
succeeded in explaining the high pion multiplicity; they suppose that annihilation 
takes place only between the cores of nucleons; the pions resulting from the annihilation 
have to be added to those coming from the virtual clouds wich have lost their centre 
to associate with. But we remark that the authors have derived the energy released 
in the core annihilation from the mean experimental energy of the pions that come 
out from the entire process, thus assuming in their phenomenological treatment an 
experimental result that they should explain. 

(1) E. FERMI: Progr. Theor. Phys., 5, 570 (1950); Phys. Rev., 92, 452 (1953). 

(2) R. E. MARSHAK: VII Conference of Rochester, chap. X (1957), p. 27. 
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Our attempt is to explain the high pion multiplicity by using Dyson’s 
hypothesis (*) introduced for justifying the second maximum in x-p scattering 
cross-section. Following DYSON we suppose the existence of a 7-r resonant 
interaction in the state of total isotopic spin 7=0 at a relative momentum of 
250 MeV/c and /=0. We must observe that this hypothesis is not too well 
defined, but we noticed, in making our calculations, that statistical weights 
for various final states were not too much sensitive to variations of the rela- 
tive momentum. 

For introducing this resonance effect in the frame of statistical theory, we 
have adapted to our problem a method suggested by LANDAU and developed 
by BELEN'KJI (6) (see Appendix). 

Numerical computations of statistical weights were made by means of 
saddle point approximation, that is exactly relativistic (7). 

In our calculations we took in account conservation of total isotopic spin, 
and assumed for interaction radius the Compton wavelength of 7-mesons. 

In Table I the statistical weights S(n_) are reported for various final 
states (*). 


TABLE I. — Statistical weights for various final states (K-meson production is not considered). 


| 
ne | S(n_) 

€ | È 

2 | 3.10% 

3 18.06% 

4 34.85%, 

5 24.29%, 

6 19.7 % (+) 

n= 4.4 

| (*) The value here reported for S(6,) is not exactly 
| calculated; this approximation does not infirm the con- 
| clusions of this work. 


(?) W. H. BARKAS, R. W. BircE, W. W. Caupp, A. G. EKSPONG, G. GOLDHABER, 
S. GOLDHABER, H. H. HECKMAN, D. H. PERKINS, J. SANDWEISS, E. SEGRÈ, F. M. SMITH, 
L. van Rossum, E. AMATDI, G. BARONI, C. CASTAGNOLI, C. FRANZINETTI and A. MAN- 
FREDINI: Phys. Rev., 105, 1037 (1957); A. G. Exspong: VII Rochester Conference, 
Chap SE AMS Te AS 

(4) I. I. POMERANGUK: Dokl. Akad. Nauk, 78, 88 (1951); E. C. G. SUDARSHAN: 
Phys. Rev., 103, 777 (1957); R. E. MARSHAK: VII Rochester Conference, chap. X 
(1957), p. 27. 

(5) F. J. Dyson: Phys. Rev., 99, 1037 (1957). 

(6) L. LANDAU: Statistical Physics, Sect. 75 (Gostehizdat, 1951); S. Z. BELEN'KJI: 
Nucl. Phys., 3, 259 (1956). 

(*) All the statistical weights reported in this work are calculated for pp (or nn) 
annihilations; the results for pn (or np) annihilations are substantially the same. 
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So we obtain the value n_= 4.4; this value is not in disagreement with 
the experimental results. 

Another experimental feature of annihilations processes is the low per- 
centage of K-mesons. In this respect we can make two hypotheses: 


a) The radius of interaction for K-mesons is given by the Compton 
wavelength of the 7-mesons. 


b) The radius of interaction for K-mesons is given by the Compton 
wavelength of the K-mesons themselves (*). 


For case a) the results are reported in Table II (+); for case b) the results 
are reported in Table III (+). 


(DABTE ALT Tics, INDI 

Ma DIE S(n_, Mx) Nw Nx S(m_, Mx) 
3 0 170% 3 0 17.93% 
4 0 | 26.05% 4 0 34.61% 
5 | 0 | 18.50% | 5 0 24.12% 
6 | 0 15 D) 6 0 19.55% (7) 
0 | 2 | 2.04% 0 2 0.06% 
Il 2 | Deis 1 2 0.18% 
D) 2 | 10 Ye (65) 2 2 028% (@) 
3 2 DAG GE (>) 3 3 ol 09) 
(*) See note(*) in Table I Percentage of cases in which a couple 
Percentage of cases in which a couple of of K-mesons appears: 0.7%. 

K-mesons appears: ~ 24%. Pion multiplicity: n= 4.4. 
Pion multiplicity: nx= 3.8%. 


Statistical weights for various final states 
Statistical weights for various final states. Interaction radius for z-mesons = #/m_e. 
Interaction radius = h/m_e. Interaction radius for K-mesons = %/mxc. 


Experimental results give (7) a value less than ~ 6% for the percentage 
of cases in which a couple of K-mesons appears. 

We can see that the case b) is the most reliable (the discrepancy between 
the value here calculated and the experimental one is not really significant 
because one can assume for the radius of interaction for K-mesons a slightly 
larger value than the K-mesons wavelength. 


(*) This hypothesis has been introduced by BARASENKOV et al. (9). 
(*) We took in account strangeness conservation, that is associated production. 
(7) G. E. A. FraLHo: Phys. Rev., 105, 328 (1957). 
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In conclusion we can say that, taking in account a rx-x resonant inter- 
action (for which there is already some evidence), the statistical theory gives 
good results for the pion multiplicity, and that the low percentage of K-mesons 
can be explained by assuming the hypothesis already introduced in the study 
of nucleon-nucleon collisions (5), namely that their volume of interaction is 
much smaller than the one of 7-mesons. 


I must expressely thank Prof. M. CixI for his continuous encouragement 
ad Dr. A. AGODI and B. VITALE for their interest in this work. 


APPENDIX 


The formula we used for evaluating the statistical weights of Table I is 
the following (5): 


S (x) =Knal Soir) + (21+1)S(n—2)x; 1Wo) + K(21+1)?S((m—4)n; 2Wo) +...]. 


Here $(n_) is the statistical weight for a final state in wich n z-mesons 
are present; So(nx); So((n — 2)z, 1W,), ..., are the statistical weights evaluated 
in absence of interaction (and not comprehensive of the coefficient that takes 
into account the conservation of total isotopic spin; this coefficient (X,,,) is 
the same for all the S,, for a given multiplicity; the various K,, have 
been calculated by YEIVIN-DE SHALIT (1)), in the cases respectively, in which 
no couple of interacting pions is present, one couple of interacting pions is 
present, and so on. 

These statistical weights have to be calculated as if in the final state 
were present: n pions, (n —2) pions and a «particle» W,, (n— 4) pions 
and two «particles » W,, and so on. 

W, has to be conceived, formally, as a particle of mass: 


where £, is the kinetic energy of the two pions at resonance in the c.m. system. 

The coefficients K,, K,,... take into account the various modes of selecting 
a couple (or two, three, ...) of interacting pions and also the fact that such 
a couple interacts only in a state of total isotopic spin 7 = 0. 


(8) A. G. Exspone: VII Rochester Conference, chap. X, p. 18. 

(°) V. S. BARASENKOV, B. M. BARBASEV, E. G. BuBELOV and V. M. MAKSIMENKO: 
Nucl. Phys., 5, 17 (1958). 

(10) Y. Yervin and A. DE SHALIT: Nuovo Cimento, 1, 1146 (1955). 
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In to order evaluate these coefficients we used those calculated by NIKISHOV (°!) 
together with obvious combinatorial methods; the results are collected in the 
following table: 


n_ | K, | Ho, 
3 | . 1.299 = 
4 2.886 1.443 
5 | 4.120 2.580 


Statistical weights reported in Table II are evaluated in a similar way. 

With regard to statistical weights of Table III, they are different from 
those of Table II only in the space weight factors of the various S,; we 
remember that the space weight factor, in the case of Table II, has the fol- 
lowing form: 


(1) 


where G is the well known factor taking into account identity of particles 
and © is the interaction volume. 

In the case considered in Table III (the interaction volume of K-mesons 
is different from that of r-mesons), (1) is replaced by (°) an average space 
weight factor: 


dae a 
n is Lee kVK AS 


n + k 


il 
= ; I k—1 yn+k-1 
ea Va ; 


where n» is the number of pions, k is the number of K-mesons, V_ and Vx 
are the space weight factors (1) for pions and K-mesons respectively, and 
E = Vx/V; = (m,/mx)*. 


(11) A. I. Nixisnov: Journ. Exp. Theor. Phys., 3, 976 (1956, American Edition). 


RIASSUNTO 


Si è calcolata la molteplicità pionica nell’annichilazione di un antinucleone ed 
un nucleone introducendo un’ipotesi proposta da Dyson riguardante una inte- 
razione di risonanza pione-pione. Si è calcolata inoltre la percentuale di mesoni K pro- 
dotti supponendo che il raggio del loro volume d’interazione sia dato dalla lunghezza 
d’onda Compton del mesone K stesso. I risultati così ottenuti sono in buon accordo 
con quelli sperimentali. 
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The Use of the Neon Flash-Tube 
for the Precise Location of Particle Trajectories. 


F. ASHTON, S. KISDNASAMY and A. W. WOLFENDALE 


Department of Physics, The Durham Colleges in the University of Durham 


(ricevuto il 18 Gennaio 1958) 


Summary. — The passage of fast cosmie ray particles through an array 
of neon flash-tubes has been studied and the accuracy of location of 
the particle trajectories has been found to be as high as in conventional 
cloud chambers. 


1. — Introduction. 


A new technique for the detection of ionizing particles has been reported 
by CoNnvVERSI et al. (*) and further developments have been made by BAR- 
SANTI et al. (2) and GARDENER et al. (*). The basic detector is a glass tube filled 
with neon which, under certain conditions, can be made to glow after traversal 
by an ionizing particle. It has been pointed out by GARDENER et al. that such 
a detector should make a satisfactory particle-locator in a magnetic spectro- 
graph. In particular, the device should be invaluable in cosmic ray studies 
where detecting layers having large areas of particle collection and high spatial 
resolution are required. 

Before such a new technique can be used in this way a detailed study of 
various aspects is necessary: this paper describes such a study. 


(1) M. Conversi and A. Gozzini: Nuovo Cimento, 2, 189 (1955). 

(2) G. BARSANTI, M. CONVERSI, S. FocARDI, G. P. Murtas, C. Ruppia and G. To- 
RELLI: Proceedings of C.E.R.N. Symposium, II, 56 (1956). 

(3) M. GARDENER, 8. Kispnasamy, E. RòssLe and A. W. WOLFENDALE: Proc. 
Phys. Soc., B70, 687 (1957). 
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2. — Experimental arrangement. 


In constructing a detecting array the first problem is the optimum tube 
diameter to choose. For tubes having 100% internal detection efficiency and 
infinitesimal wall thickness, the number of tubes required to cover a given 
area and produce a given spatial resolution is virtually constant, viz: 6 layers 
of tubes of diameter d give nearly the same resolution as 2 layers of tubes 
of diameter d/3. However, considerations of the gas pressure required to give 
a reasonable number of ion-pairs and the difficulties involved in manipulating 
tubes having very thin walls dictate a tube diameter of several mm. The 
tubes used in the experiments to be described had internal diameter 5.9 mm 
and external diameter 7.7 mm. 

The apparatus used to study the accuracy of location of particle trajectories 
is shown in Fig. 1. Four tube arrays, A, B, C and D, and associated Geiger 
counters served to trace the trajectories of high energy particles at sea level. 
The tubes were rigidly mounted 8.0 mm apart between parallel electrodes: 
scattering was reduced by replacing the central regions of the electrodes by 
Aluminium foil. The tubes were filled with commercial neon at a pressure 
of 65 cm Hg and the residual air pressure was estimated to be 10-? mm Hg. 
The characteristics of the high voltage pulse applied across the plates after 
a 4-fold coincidence were as follows: peak height 6.3 kV/cm, time delay to 
the start of the pulse 2.0 us, rise-time to 75% of the peak height 0.5 us, pulse 
length 4.0 us. Under these conditions the internal efficiency of the tubes was 
70% and the layer efficiency 52%. In order to reduce the scattering of particles 
in the layers a low momentum cut-off was imposed by demanding further 
penetration of a lead block below the apparatus. When operated as a spectro- 
graph a magnet will be inserted between arrays B and C and the incident 
and emergent directions determined by AB and CD respectively. In the expe- 
riment reported here the magnet was removed but the angle, 0, between the 
apparent trajectories in AB and OD was measured; the distribution of this 
angle after correction for scattering gave the noise-level or « no-field » limit 
of measurement. 

Before proceeding it is necessary to consider the optimum geometrical 
arrangement of the individual layers of tubes within an array. Tf the tubes 
are arranged vertically one above the other it is possible for vertical particles 
to pass through the dead space between the sensitive regions of adjacent tubes 
in each layer with consequent loss in overall efficiency of detection, since 
events where no tube flashes in one array cannot be used. The overall efficiency 
is defined as the ratio of the number of measurable photographs to the total 
number taken. Similarly a vertical particle passing through the sensitive 
region of a tube in each layer can give rise to considerable uncertainty in the 
location of the trajectory. A staggering angle can be chosen to minimize these 
effects if the angular distribution of the particles is known and this has been 
done here. The practical performance of this arrangement has been compared 
with that for vertical arrays and whereas the overall efficiency of particle 
detection is higher the increase in precision of location in the former case is 
not very great. 
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Fig. 1. — The experimental arrangement. (a) Shows the arrangement of 4 tube arrays 

(A, B, C, D) and the cotton fibres used to align the arrays (A, B and C, D are photo- 

graphed by separate cameras). (b) Shows the passage of a single fast y-meson through 
the apparatus. 


5 40 — Il Nuovo Cimento. 
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3. — Measurements of the trajectories. 


Measurements of the track directions in AB and CD were made by 3 methods: 


(i) The centre of gravity of the array formed by the centres of the tubes 
which flashed at each level was calculated and the inclination of the line joining 
the appropriate pairs was determined. 


(ii) A full-seale diagram of the tube assembly was drawn and a cotton 
fibre adjusted over its surface to give the best estimate of the track direction. 


(iii) The photographic records were projected on to a rotable screen 
ruled with close parallel lines and the direction of a line determined which 
satisfied the condition that the sum of the distances from the centres of the 
flash-tube images to the line was a minimum. 


In each method where a tube was observed to flash in such a position as 
to be incompatible with the trajectory defined by the other tubes this flash 
was rejected. Several comments can be made on these methods. Method (i) 
has the advantage of being completely objective but some information present 
in the record is not used, that is the reduction in uncertainty of position at 
one level arising from the knowledge of the approximate direction given by 
the information from the two levels taken together. In principle this can be 
overcome by applying the method of least squares, say, to the data from 
the two levels in turn, but the calculations are tedious. 

Method (ii) is subjective but has the advantage that not only the excited 
tubes, but also the gaps between the sensitive regions, can be used to locate 
the trajectories. 

Method (iii) on the other hand relies only on the excited tubes since these 
are the only images visible on the film. A basic assumption is that the image 
accurately represents the end of the tube, and with carefully constructed tubes, 
this condition is attained. 

All the photographs taken in this experiment have been measured using 
methods (ii) and (iii) and a smaller number with method (i). The conclusion 
is that method (i) is inferior to the other two and that there is little difference 
between the results of methods (ii) and (iii); method (iii) is quicker to apply 
and is indicated for applications where a high particle rate is encountered. 
Where great significance is attached to individual particles, checking by 
method (ii) is desirable. A further check can be carried out by determining 
the points of intersection with the central plane of the trajectories found from 
the upper and lower halves of the apparatus. The separation of these points 
will be small for single particles. 

The distribution in 6 found from measurements with the present apparatus 


has been found to be 0.33°+0.02°. This distribution arises from the fol- 
lowing causes: 


a) The scattering of particles in arrays B and 0. 


b) The error of setting the measuring device (cross-wire etc.) on the 
best estimate of the trajectory. 
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c) The error of interpretation, i.e. the difference between the estimated 
and true trajectories. 


d) The traversal of the apparatus by more than one particle. 


Of these causes only c) and d) need further amplification. Item c) refers 
to cases where one or more of the tubes discharged were not traversed by the 
ionizing particle itself. Probably the most important contribution arises from 
unrecognized knock-on electrons, e.g. a case where the fast particle passes 
through the insensitive region of one tube and an associated knock-on electron 
passes through the next tube such that a systematic difference between the 
best estimate and the true trajectory cannot be recognised. Also included 
in c) is the spread of true trajectories about those best-estimates derived using 
the criterion of minimum sum of the distances from estimated trajectory to 
the centres of the excited tubes. 

The contribution from scattering and errors of setting are known and thus 
the error from the combination of e) and d) can be determined. Writing o.., 
Os, Orm for the r.m.s. deflections due to scattering, setting, and the remaining 
causes, we have 


Ou = 0.22, of —0.14° giving Gen =0.20°+0.03° 


It is usual to quote the uncertainty as that of position at a single level (e.g. the 
middle of each array) and this can be calculated. It must be borne in mind 
however that this uncertainty is less than would be observed if the single 
level alone were observed on account of the information on direction given by 
the combined data from the two levels. 

The required r.m.s. uncertainty, 0, is given by 6 = Loyem/2 where L is the 
Separation of each pair of measuring levels. The result is 


Ô = (0.62+0.10) mm. 


This value can in principle be compared with that expected from the geo- 
metry of the arrangement and the efficiency characteristics of single tubes. 
Limiting forms of the efficiency variation across the tube diameter have been 
considered and the value of 6 determined graphically in each case. The result is 


(6; = 0.46 mm, 


for an efficiency function 7 = f(r/a) where a is the radius and r is the distance 
from the centre, with f(r/a) = 1.0 for r/a < 0.70 and f(r/a)= 0 for r/a > 0.70 and 


for f(r/a) = 0.70 for all values of r/a<1 

The actual form of the efficiency variation across the tubes was studied 
in a separate experiment in which the arrangement of arrays was modified. 
Array A was placed mid-way between arrays B and C and tracks were defined 
by the observed flashes in B and C. The distances of these tracks from the 
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contres of those tubes in A through which the particles had passed were 
accurately measured. The resulting variation of internal efficiency with di- 
stance from the centre is shown in Fig. 2. Also shown is the expected variation 
according to a theoretical analysis of LLOYD (to be published) where nor- 
malization to the same internal efficiency has been made. A broadening of 
the expected distribution is necessary due to the uncertainty in track location 
(— 0.3 mm) but this is too small to be significant. It is apparent that the 
agreement with theory is good. 

Since the form of the efficiency vari- 


ino = ation is roughly mid-way between the two 
de Theoretical extreme assumptions mentioned above 
80 È > PILA it is reasonable to expect that the re- 
sulting uncertainty be mid-way between 
60 ; 6, and 6,. The observed value, of 
0.64 mm, is thus in good agreement with 

20 expectation. 

20 

Fig. 2. — The variation of internal effi- 


ciency with distance of the ionizing particle 

from the centre of the tube. 7 is the inter- 

nal efficiency and r the distance from the 

tube centre. The theoretical variation allows for the loss of electrons by diffusion to 

the walls in the time delay occurring between the passage of the particle and the 
application of the high voltage pulse. 


4. — Conclusions. 


It can be concluded from the discussion of the accuracy of location that 
the effect of spurious phenomena such as discharges produced by unrecognized 
knock-on electrons, is small. 

A comparison can be made between this method and the cloud chamber. 
In a cloud chamber the inherent track width, for cosmic ray studies, is ~ 1 mm. 
It is possible to determine the centre of a track to within about 0.2 mm but 
to this uncertainty must be added that arising from gaseous distortion. By 
careful temperature control distortion can be reduced but both long and 
short term variations in distortion produce difficulties. In the Manchester 
high energy spectrograph (HOLMES et al. (*)) cloud chambers were used as 
detectors and the overall uncertainty in track location was shown, in a sub- 
sidiary experiment by LLOYD et al. (5) to be about 0.8 mm. Further reduction 
by more careful temperature control is presumably possible but the location 
uncertainty in a cloud chamber is inherently variable and in that respect 
alone is inferior to the flash-tube array with its constant geometrical uncertainty. 


(4) J. E. R. HoLmrs, B. G. Owen, A. L. RopGERSs and J. G. WiLson: Proc. Phys. 
Soc., A 68, 793 (1955). 


(5) J. L. Luoyp, E. RòssLe and A. W. WoLFENDALE: Proc. Phys. Soc., ATO, 
421 (1957). 
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It is apparent that a further reduction in error of location can be made 
in a flash-tube array by increasing the number of layers, an upper limit being 
set by the maximum tolerable scattering uncertainty. In cosmic ray expe- 
riments, however, the limit is liable to be set more by the lack of a sufficient 
flux of particles of small magnetic deflection than the difficulty in measuring 
the deflections. 


The authors wish to thank Professor G. D. RocHESTER for his interest 
in the work and for useful suggestions. We are grateful to our colleagues 
Messrs. M. GARDENER, D. G. JoNES and J. L. LLOYD for their assistance in 
many ways. 

The initial development work on the detectors was supported by the 
Department of Scientific and Industrial Research and two of us (F.A. and S.K.) 
wish to thank this body for the provision of a Research Studentship and 
Research Assistantship respectively. The apparatus itself was constructed 
with a grant from the Royal Society. 


TILA SING) (5) 


E stato studiato il passaggio di radiazione cosmica veloce attraverso un rivelatore 
costruito con scintillatori al neon e l’esattezza dell’individuazione delle traiettorie delle 


particelle è risultata altrettanto buona di quella permessa dalle camere a nebbia con- 
venzionali. 


(*) Traduzione a cura della Redazione, 
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Study of +'-Meson Decay in a Propane Bubble Chamber (*). 
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University of Michigan - Ann Arbor, Michigan 
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University of Michigan - Ann Arbor, Michigan 
Brookhaven National Laboratory - Upton, New York 


(ricevuto il 22 Marzo 1958) 


Previous large sample studies of 
t+-meson decay events have been carried 
out by means of the nuclear emulsion 
technique exclusively (1). There is, of 
course, always the possibility of a geo- 
metric bias existing in experiments of 
this sort and therefore it is of interest 
to examine data which have been col- 
lected and analyzed by a different expe- 
rimental method. 

In a recent exposure of the Uni- 
versity of Michigan propane bubble 
chamber to a K+-meson beam at the 
Brookhaven Cosmotron, 136 7*-mesons 
were observed to decay from rest. The 
details of the method of collection and 
processing of the data have been con- 
sidered elsewhere and will not be dis- 
cussed here (2). The purpose of this note 


(*) Supported in part by the U.S. Atomic 
Energy Commission. 

(*) Present address: U.S. Naval Research 
Lab., Washington 25, D.C. 

(1) Proceedings of the Seventh Annual Ro- 
chester Conference, 7, 19 (1957). 

(ND SC MRYER, Wily Ibi, Tidy hah D AN, 
GLASER: Phys. Rev., 107, 279 (1957). 


is to present, briefly, the results of the 
analysis of the t+ decay events. 

For each event the ranges of the 
stopping pions were measured. In ad- 
dition, the angles included between the 
pion tracks were also obtained. This 
information, of course, overdetermines 
the event kinematically, but does have 
the advantage that it allows a check on 
the internal consistency of the experi- 
mental method. As a further check 
over a hundred u+-meson tracks arising 
from 7’s decaying from rest were also 
measured. The average u+ length ob- 
tained from these measurements was 
found to be in good agreement with the 
range predicted by the range energy 
tables calculated for propane (C3H,). 

For 118 of the decays the energies 
and the charge of the three pions was 
unambiguously determined. The re- 
maining 18 events were such that the 
T -meson and one of the z+-mesons 
escaped from the chamber. Thus, in 
these cases the charge of the two 
escaping particles is ambiguous. It is 
to be noted that it is only the charge 
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TaBLe I. — The logarithms to the base 10 of the probability ratios are given for the 


Les 


and 2+ spin-parity combinations compared to the 0° spin-parity combination, 


as calculated by the Dalitz-Fabri analysis. For the expected values the standard deviat- 
ions are given. 


log. Pol flos;g21a/P0x | los, baa) 0 

Observed value if 18 ambiguous 

events are distributed as 07 . — 13 — 28 — 11 
Observed value if 18 ambiguous 

events are distributed as other 

spin-parity combination. . — 12 — 23 — 8 
Expected value if 07 is correct — 1204 — 539 + 6 — 49 + 13 
Expected value if other spin- 

parity combination is correct + 8+2 + 133 + 30 +39 + 8 


which is ambiguous since the energies 
of the escaping particles are known from 
the angular measurements. 

The energy distribution, Fig. 1, and 
the cos 0 distribution, Fig. 2, are the 
same within statistics as the distrib- 
utions derived from emulsion measure- 
ments (1). Therefore no systematic bias 
in the emulsion method has been dis- 
closed by the bubble chamber method. 
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Fig. 1. — Distribution of the events versus the 

energy of the 7 -meson. e is the kinetic energy 

of the x -meson divided by the maximum kinetic 
energy of the = -meson. 


The Dalitz-Fabri analysis of the data 
supports the 07 spin-parity combination, 
as can be seen from Table I which 
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presents the probability ratios. The 
agreement of the probability ratios with 
the 0 spin-parity combination cannot 
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Fig. 2. — Distribution of the events versus cos 0. 

8 is the angle in the center of mass system 

between the momentum vector of the 77 -meson 

and the relative momentum vector of the two 
n*-mesons. 


be altered by distributing the 18 am- 
biguous events in the most unfavorable 
way. A y? test on the energy spectrum 
yields a value of P = 25% when com- 
pared to the 0° distribution function 
and for the angular distribution a value 
of P= 7% is obtained. A complete list 
of the events can be obtained from the 
first author (T.F.Z.). 


LIBRI RICEVUTI E RECENSIONI 


B. E. DuckwortH - Mass spectro- 
scopy. Cambridge Univ. Press., 
1958, pp. XVI - 206. 


Il libro, che fa parte dei « Cambridge 
Monographs on Physics », riassume in poco 
meno di 200 pagine l’attuale situazione 
nel campo della spettrometria di massa 
e delle sue principali applicazioni. 

L'esposizione del DucKwor tH, chiara 
ed ordinata, illustra concisamente 1 pro- 
gressi compiuti e i dati acquisiti. I pro- 
blemi fisici sono trattati esaurientemente, 
ma senza entrare nei dettagli, in modo 
da rendere il libro di facile lettura al non 
specialista; in questo stesso spirito le 
questioni più tecniche, come l’elettronica 
e l’alto vuoto, sono state evitate per non 
appesantire il libro. 

Dopo una introduzione storica, cinque 
capitoli del libro sono dedicati all’ottica 
ionica e alla strumentazione, due alla de- 
terminazione rispettivamente delle masse 
atomiche e delle abbondanze isotopiche 
e tre alle applicazioni. Per quanto ri- 
guarda queste ultime, solo quelle riguar- 
danti la fisica e la geologia sono trattate 


in qualche dettaglio: un intero capitolo 
concerne le applicazioni in fisica nucleare, 
e relativamente ampio spazio è dedicato 
alla ionizzazione e alla dissociazione delle 
molecole soggetto ad urto da parte di 
elettroni, e allo studio degli ioni meta- 
stabili; infine sono trattate le applica- 
zioni alla geologia, sia quelle riguardanti 
gli isotopi radiogenici (determinazione 
delle età geologiche con i metodi U-Pb 
A-K e Rb-Sr), sia quelle riguardanti gli 
isotopi stabili (variazioni delle abbon- 
danze in conseguenza di processi chimico- 
fisici occorrenti in natura). Vengono solo 
sfiorati i problemi in cui lo spettrometro 
di massa funziona esclusivamente come 
strumento analitico, per esempio l’uso 
degli isotopi per lo studio del mecca- 
nismo delle reazioni chimiche e i metodi 
di separazione degli isotopi. 

Il libro è di qualche interesse anche 
per lo specialista; in particolare la biblio- 
grafia è molto ampia e, per quanto ri- 
guarda la strumentazione, pressochè com- 
pleta. 
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